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(III) Status of Claims 

Claims 1-10 and 12 are pending and were all finally rejected in an Office Action dated 
May 4, 2004 ("the Action"). Claims 1-10 and 12 are the subject of the present appeal. A copy 
of the appealed claims is attached to this brief in a Claims Appendix. 

(IV) Status of Amendments 

No amendments have been filed subsequent to the final rejection. 

(V) Summary of Claimed Subject Matter 

Independent claim 1 is directed to a method for identifying at least one binder molecule 
that interacts with a polypeptide or protein. The claimed method recites six steps, (a) through (f) 
(specification at page 2, line 26, to page 3, line 16). In step (a), two amino acid types (AAl and 
AA2) in a polypeptide or protein of interest are chosen, whereby AA2 at least once occurs 
directly subsequent to AAl in the amino acid sequence of the polypeptide or protein 
(specification at page 2, line 28, to page 3, line 2). In step (b), the two amino acid types in the 
polypeptide or protein are labeled with ^^C (AAl) and ^^N (AA2) (specification at page 3, 
lines 3-5). In step (c), a first NMR spectrum of the labeled polypeptide or protein is generated, 
which spectrum identifies signals fi"om the labeled amino acid pair (specification at page 3, 
lines 6-8). In step (d), the labeled polypeptide or protein is contacted with a potential binder 
molecule or a mixture of potential binder molecules having a molecular mass of firom 50 to 
1000 Da (specification at page 6, lines 19-21). Contacting step (d) is carried out under 
conditions and for sufficient time so as to allow binding between the polypeptide or protein and 
the potential binder molecule(s) (specification at page 3, lines 9-1 1). In step (e), an NMR 
spectrum fi-om the mixture of step (d) is generated (specification at page 3, lines 12-13). In 
step (f), the first and second NMR spectra are compared to identify a chemical shift change 
between the two spectra, wherein a shift indicates an interaction between the potential binder 
molecule and the labeled polypeptide or protein (specification at page 3, lines 14-16). 

Dependent claim 12 fiirther limits claim 1 by requiring that the method be used for 
screening a compound library (specification at page 4, lines 12-13). 
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(VI) Grounds of Rejection to be Reviewed on Appeal 

(A) Claims 1, 4, 5, 8-10, and 12 stand rejected under 35 U.S.C. § 103(a) as allegedly 
unpatentable over Yabuki et al. (1998) J. BIOMOLECULAR NMR, 1 1 :295-306 (hereinafter 
"Yabuki") in view of Moore, U.S. Patent No. 6,060,603. 

(B) Claims 2, 3, 6, and 7 stand rejected under 35 U.S.C. § 103(a) as allegedly 
unpatentable over Yabuki in view of Moore and in further view of Fesik et al., WO 97/18471 
(hereinafter "Fesik"). 

(VII) Argument 

(A) Rejection of claims K 4, 5, 8-10, and 12 under 35 U.S.C. S 103fa) as allegedlv unpatentable 
over Yabuki in view of Moore 

To establish a prima facie case of obviousness under 35 U.S.C. § 103(a), there must be 
some suggestion or motivation, either in the cited references or in the knowledge generally 
available to one of ordinary skill in the art, to modify a single reference or to combine it with one 
or more other references. The suggestion or motivation to make the claimed combination must 
be found in the prior art and not based on applicant's disclosure. In re Vaeck, 947 F.2d 488, 
20 USPQ2d 1438 (Fed. Cir. 1991); MPEP § 706.02(j). 

Yabuki describes a dual amino acid-selective labeling method and its use for 
investigating structures of proteins and protein complexes. Yabuki does not suggest using the 
labeling technique to carry out a screen to identify whether a candidate molecule has the ability 
to bind to a given target protein. Instead, Yabuki describes the use of the labeling technique to 
investigate the nature of the interaction between two proteins (Ras and Raf) that were previously 
known to bind to each other. 

The Action concedes that Yabuki does not disclose carrying out its labeling method with 
a binder molecule having a molecular mass of from 50 to 1000 Da (the Action at page 3). 

Moore describes methods for modeling biologically active ligands and designing 
mimetics of such ligands. The Action alleges that it would have been obvious to one of ordinary 
skill in the art to modify the disclosure of Yabuki to include the use of small molecular weight 
ligands taught by Moore (the Action at page 4). 
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The person of ordinary skill in the art, at the time the present apphcation was filed, would 
have lacked the requisite suggestion or motivation to modify Yabuki in the manner contemplated 
by the Action. Moore does not describe any molecule having a molecular mass of from 50 to 
1000 Da that has the ability to bind to Ras or any other protein described in Yabuki. 
Furthermore, nothing in Moore suggests modifying the protein structure analysis method of 
Yabuki to instead carry out a screening method that compares a first NMR spectra (of a labeled 
polypeptide or protein) and a second NMR spectra (of the labeled polypeptide or protein that has 
been contacted with a potential binder molecule having a molecular mass of from 50 to 
1000 Da). 

Because Yabuki is concerned with investigations of protein structure and the 
characterization of protein complexes, the skilled artisan would have had no reason to modify the 
methods of Yabuki to create a method that screens candidate compounds to evaluate their ability 
to bind to Ras or any other protein. Furthermore, neither Yabuki nor Moore describe a single 
low molecular weight compound (i.e., having a molecular mass of from 50 to 1000 Da, as is 
required by the claims) that was known to bind to Ras that the skilled artisan even might have 
attempted to use in the method of Yabuki in place of Raf to investigate the nature of a Ras-ligand 
interaction. 

The Action stated that "absent evidence to the contrary, the range recited in the instant 
claims from 50-1000 Daltons is viewed as mere optimization of the prior art assay" (the Action 
at page 4). Because Yabuki does not describe a screening assay to identify potential binder 
molecules (but instead describes experiments that evaluate Ras binding to Raf), there would have 
been no basis or rationale for the skilled artisan to "optimize" the methods of Yabuki by, for 
example, replacing Raf (a known hgand of Ras) with some other molecule having a molecular 
mass of from 50 to 1000 Da. Yabuki describes a dual amino acid-selective labeling method 
using a protein (Ras) and a known binding partner (Raf) and thus gives no suggestion or 
motivation to replace Raf with low molecular weight compounds having no known Ras-binding 
ability. 

Accordingly, Appellants respectfully submit that the Action does not establish a 
prima facie showing of obviousness and request that the rejection of independent claim 1 and 
dependent claims 4, 5, 8-10, and 12 be withdrawn. 
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(B) Rejection of claims 2, 3. 6. and 7 under 35 U,S.C. § 103(a) as allegedly unpatentable over 
Yabuki in view of Moore and in further view of Fesik 

The Action contends that it would have been obvious to one of ordinary skill in the art to 
incorporate a comparison method of the various assays as taught by Fesik into the method of 
Yabuki in view of Moore to compare the binding of ligands to various biomolecules determined 
by NMR and to also observe chemical shifts detected by 2-D NMR techniques (the Action at 
page 5). 

As detailed in Section (A) above, Yabuki describes the use of a dual amino acid-selective 
labeling method to investigate structures of proteins and protein complexes. Nothing in Yabuki 
suggests generating the two NMR spectra recited in the claims and subsequently comparing the 
two spectra to identify an interaction between a potential binder molecule having a molecular 
mass of from 50 to 1000 Da and a labeled polypeptide. As detailed above, Moore provides no 
suggestion or motivation to modify the technique of Yabuki to convert the method described 
therein into an assay to screen for interactions between a labeled polypeptide and a potential 
binder molecule having a molecular mass of from 50 to 1000 Da. Accordingly, the combination 
of Yabuki and Moore do not render obvious claim 1 or the claims that depend therefrom. 

Fesik does not add what is lacking in Yabuki and Moore. Fesik' s disclosure of processes 
for identifying compounds would not have provided the requisite suggestion or motivation to 
cause the skilled artisan to so thoroughly alter the methods described by Yabuki so as to, instead 
of characterizing protein structures as taught in Yabuki, modify the labeling techniques described 
therein to be used for the altogether different purpose of carrying out a screening assay to 
identify potential binder molecules. Yabuki is not directed to screening assays and nothing in 
Moore or Fesik, considered either alone or in combination, would have led the skilled artisan to 
(a) carry out the methods of Yabuki with a molecule having a molecular mass of from 50 to 
1000 Da, and/or (b) fundamentally modify the techniques of Yabuki so as to screen to identify 
binder molecules. 

In light of these comments, Appellants respectfully submit that the Action does not 
establish a prima facie showing of obviousness and request that the rejection of dependent 
claims 2, 3, 6, and 7 be withdrawn. 
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(O Rejection of claim 12 under 35 U.S.C. $ 103(a) as allegedly unpatentable over Yabuki in 
view of Moore 

Claim 12 was not rejected under a separate heading in the Action. However, appellants 
consider claim 12 to be patentable for the following separate reasons, in addition to those reasons 
provided in Section (A) above. 

Dependent claim 12 requires that the method of independent claim 1 be used for 
screening a "compound library." Nothing in Yabuki or Moore, taken alone or in combination, 
suggests applying the dual amino acid-selective labeling method of Yabuki to a "compound 
library" to identify molecules in the library that bind to a polypeptide or protein of interest. Li 
particular, there is no description whatsoever in Moore of any compound library that the skilled 
artisan would have had reason to use in a labeling method of Yabuki. 

Accordingly, Appellants respectfully submit that the Action does not establish a 
prima facie showing of obviousness and request that the rejection of dependent claim 12 be 
withdrawn. 

For the foregoing reasons. Appellants respectfully submit that the Action's rejections of 
claims 1-10 and 12 are unwarranted. Appellants therefore request that the Board overturn the 
rejections. 

Respectfully submitted. 

Date: ^^7,iOOV^ 

Fish & Richardson P.C. 
Citigroup Center 
52nd Floor 
153 East 53rd Street 
New York, New York 10022-461 1 
Telephone: (212) 765-5070 
Facsimile: (212) 258-2291 

30211243.doc 




Reg. No. 47,443 
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Claims Appendix 



1. A method for identifying at least one binder molecule comprising the steps of: 

(a) choosing two amino acid types (AAl and AA2) in a polypeptide or protein of 
interest, whereby AA2 at least once occurs directly subsequent to AAl in the amino acid 
sequence of the polypeptide or protein, defining an amino acid pair AA1-AA2; 

(b) labeling the two amino acid types (AAl and AA2) in the polypeptide or protein of 
interest, whereby all AAl -residues is labeled with ^"^C and all AA2-residues with ^^N; 

(c) generating a first HNCO-type NMR spectrum of the labeled polypeptide or protein 
from step (b), thereby identifying signals from the labeled amino acid pair AA1-AA2; 

(d) contacting the labeled polypeptide or protein with a potential binder molecule or a 
mixture of potential binder molecules under conditions and sufficient time for allowing binding 
of the potential binder molecule or the potential binder molecules and the labeled polypeptide or 
protein, wherein the potential binder molecule or potential binder molecules have a molecular 
mass of from 50 to 1000 Da; 

(e) generating a second HNCO-type NMR spectrum, or a ^H-*^N correlation type NMR 
spectrum, of the mix from step (d), monitoring signals identified in step (c); and 

(f) comparing the first and the second NMR spectra, whereby a chemical shift change of 
the signals identified in step (c) between the two spectra indicates an interaction between the 
potential binder molecule or the potential binder molecules and the labeled polypeptide or 
protein. 

2. The method of claim 1, wherein the labeled amino acid pair AA1-AA2 is unique 
within a sphere radius of 10 A within the polypeptide or protein, wherein the 10 A radius is 
measured from the labeled amino acid pair AA1-AA2. 

3. The method of claim 1, wherein the labeled amino acid pair AA1-AA2 is unique 
within a sphere radius of 50 A within the polypeptide or protein, wherein the 50 A radius is 
measured from the labeled amino acid pair AA1-AA2. 
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4. The method of claim 1, wherein the labeled amino acid pair AA1-AA2 is unique 
within the polypeptide or protein. 

5. The method of claim 1, wherein the labeled amino acid pair AA1-AA2 is within a 
binding pocket of the polypeptide or protein. 

6. The method of claim 1, wherein the labeled amino acid pair AA1-AA2 is in the 
proximity of an active site within the polypeptide or protein. 

7. The method of claim 1, wherein the result of the method is compared to the result of 
any other suitable binding or activity assay. 

8. The method of claim 1, wherein the polypeptide or protein has a size of from 10 to 
150 kDa. 

9. The method of claim 1, wherein the potential binder molecule is a peptide, 
polypeptide, protein, antibody, nucleic acid molecule, or carbohydrate. 

10. The method of claim 1, wherein the potential binder molecule comprises a complex 
of one or more of a peptide, polypeptide, protein, antibody, nucleic acid molecule, or 
carbohydrate. 

12. The method of claim 1, wherein the method is used for screening a compound 

library. 



8 



Applicant 
Serial No. 
Filed 
Page 



Johan Weigelt et al. 
09/986,240 
October 19, 2001 
9 of 10 



Attorney's Docket No.: 13425-047001 / 00357-US 



Evidence Appendix 



Enclosed are copies of the following evidence relied upon by the Examiner as to grounds 
of rejection to be reviewed on appeal. 

Exhibit 1 : Yabuki et al. (1998) J. BIOMOLECULAR NMR, 1 1 :295-306 
Exhibit 2: Moore, U.S. Patent No. 6,060,603 
Exhibit 3: Fesik et al, WO 97/18471 
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Related Proceedings Appendix 

There are no appeals or interferences for related cases. 
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Dual amino add-selective and site-directed stable-isotope labeling of tbe 
human c-Ha-Ras protein by ceU-free synthesis 
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Abstract 

We developed nro »i«hods for stable-iswope labelmg of pmeiitt by ceD-fiee synihcsU. Risdy. we applied adl-ftee syndesis to the dual 
amioo aeiiselective "C-"N tabeliog mediod, oripnaUy devel<ved fac lnvim syaenu by Kafaodio and po^wwtes. For to purpose, 
we IOC* one of the advanagej of a ccD-fiw proidn synihesis system; die anuno add^iecdve sablc-jsotope labdmg k we <« ao»pe 
iSii poble«!TbeX« of «leedveXetv.don w« Th.35 seP' in d« -effecwr 3^> f^J^J^- 

comple«d Xid. d>e Ras-Wnding domain of rfUf-l (Raf RBD) (the totd mol«»lar mass is about 30 U^. tftuig a 15inL ftrtenetei cob 
edl-ftec system, which was optimized to produce about 0.4 mg of Ras protein per l-mL reaction, vwih 2 eaA rfDW Clprotos and 
L-[>5N]tbreonine, wc obtained about 6 of Ras protein. A» die Pro-Tlir sequence is unique in die protein. Ae Hif" cro» petk of *c 
R«s.Raf RBD complex was unamHeuoudy identified by the 2D >H.>5n HNCO cstpaiment. The Ser^* ooss peak >«as nmilady ideimfied 
widi dM ["C]Asprt"N]Ser-«lecih<ely labeled Ras protein. Tbete were no Isotope sciamblfag problems in diB sniiftr. Secondly, we have 
established a meduxJ for producing * uuUignnn quanriiy of siie-^pedfically stable^sotope .labeled protein by a cell«ce sysum mvolv^ 
amber suppression. Tbe £ coB amber suppicssor tRNA^^^ (25 mg) was piepaied hy inwro transcripnon widi T7 polymeiast We 
aminoecylated die OOfA^cUA transcript widi purified £ ecM tfrnyl-OtHA syndietase. osing 2 nje of L-r*»N]QTOsine. In flic gene enco<^ 
the Ras protein, d« eodon for was changed to an amber dXton CTaO). This t«iplale «nd AeTl^NlTyr^T^^^ 
Ibr 30 minis 30 mL of SeaKcdl-fieesystent The subsequent puriiieadonyieldeJ^2^n^ H- r*H5QC 

specmim of d»e labeled Ras proteia only one cross peak was observed, wWch was unambiEUouslyaisi^ . 
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IntrodactioKi 

Protein biosynthesis systems, such as the Escherichia 
coli coupled transcnption-translation system, have 
been remaxkably improved, and a milligiam quamicy 
of protein can be prepared (Kigawa ec al^ 1 995; Kam 
et aL, 1996). For NMR spectroscopy, cell-free pro-- 
tdn synthesis has a number of advantages over the 
ordinary recombinant protein synthesis in vivo. We 
have been applying the cell-free protein synthesis to 
stable-isotope labeling (Kigawa et aL, 1995). 

Kainosho and co-workers have developed a dual 
amino acid-selective ^^C-^^N labeling technique 
(Kaino^o and Tsuji, 19S2; Kainosho et aL, 1985, 
1987; Wesiler et al,, 1988a,b). This technique uti- 
lizes protein samples in wfaich the main chain caxbonyl 
caibons of one amino acid type are labeled widi '^C 
and the amide nitrogens of another amino acid type 
are labeled with ^^N. The NMR signals of the amino 
acid residues that possess a ^^CO-^^N linkage can be 
extracted on the basis of the ^^C-^^N spin coupling. 
This technique was apphed to an investigation of the 
structure of Streptomyces subtilisin inhibitor and its 
complex with the target enzyme^ subdlisin (Kainosho 
and Tsuji, 1982; Kainosho et aL, 1985. 1987; West- 
1^ et al., 1988a,b), and was also applied to structural 
analyses of an antibody and its fragment (for a review^ 
see Arau et aL, 1994). As mdicated in these stud- 
ies, the advantage of diis dual labeling technique is 
that sequence specific assignments can be performed, 
even for very large proteins, such as an ^tire IgG 
molecule (molecular mass about 150 kDa, Kato et aL, 
1989). However, if in vivo expression methods are 
used for this kind of amino acid^seleative labeling, 
then the amino acid metabolism can cause an isotope- 
scrambling problem in some cases, such as for Asp 
or Ser residues (Kainosho and Tsuji, 1982; Kainosho 
et al.. 1987; Mcintosh and Dahlquist. 1990). This is a 
serious drawback of this technique. 

This scrambfin^ problem can be solved by using 
a cell-free protem synthesis system. With this system, 
the Asp residues could be labeled without labeling the 
Asn residues, and similarly, the Ser residues could 
be labeled without labeling the Giy residues (Kigawa 
et al., 1995), whereas in vivo expression methods do 
not allow this selectivity (Mcintosh and Dahlquisi, 
1990; Yamasaki et al., 1992). In this smdy, we used 
the cell-free system for dual labeling to investigate the 
structure of the effector region of the human c-Ha-Ras 
protein complexed with the Ras-binding domain of the 
rat c-Raf-1 protein (Raf RBD; the total molecular mass 



of die complex is about 30 kDa). Tins allowed some 
of the main<hain amide signals of the e£fecior region 
to be assigned in a sequence specific manner by the 
HNCO experiments. 

In the next st^, we tried site-directed incorpo- 
ration of a stable-isotope labeled amino acid into a 
protein, which would be useful for NMR analyses. 
The site-directed incorpomdon of unnatural amino 
acids was achieved by cell-free protein synthesis in- 
volving amber suppression (Noren et al., 1989; Bain 
et al., 1991). Elhnan et aL (1992) produced T4 
lysozyme containing a unique ^^C-labeled Ala using 
chemically amincjaceylated amber suppressor OINA, 
and recorded ^^C^-filtered proton NMR spectra. The 
yield was 0. 1 5 mg of lysozyme from a 40-mL reaction 
mixoue. 

Sonar et al. (1994) and Liu ^ al. (1995) reported 
the site-directed labeling of bacteriorhodopsin wlih 
and ^^C, respecdvely; an overproduced amber sup- 
pressor, tRNA^^cUA- from£ coU was aminoacylated 
with a stable-isotope labeled L-^sine by an K coU 
crude protein fraction, and was used in ^e cell-free 
translation of a messenger RNA with an amber codon 
(UAC) at a specified site. The yield was 0.06 mg of 
bactmorhodopsin from a 2-mL reaction mixmre. 

In die present smdy, an E coli cell-free sys- 
tem was applied, to the site-directed stable-isotope 
labeling of the Ras proton. An E. coli si^pressor 
tRNA^cUA was I prepared by in vitro transcription 
with T7 RNA pdlymerase, and was aminoacylated 
with purified Qncofsyl-tRNA syndtetase. The aminoa- 
cylated cRNA^ctJA ^ tzsed in the cell-free coupled 
transcripuon-translacion of a template DNA with an 
aciber codon (TAt5) for Tyr^^. We obtained 2.2 mg 
of purified [^^NJiyr^^-Ras protein and succeeded in 
measuring the ^H-^^NHSQC spectrum. 



Materials and methods 

Preparation of the seleoiuely ^^C/^N^labeled Ras 
protein 

The Ras protein used in diis study consisted of 171 
amino add residues, and lacked the C*cerminal 18 
amino acid residues (Ha et al., 1989). The £. coli 
S30 extract for the cell-free protein synthesis was pre- 
pared according to Pratt (1984) &om£. coli strain A19 
(metB, ma). The T7 RNA polymerase was prepared 
according to Zawadzki and Gross (1991). The large 
scale reaction for selective ^^C/^^N-labeUng of die 
Ras protein was performed by our recently improved 
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protocol (to be publisbed elsewhere). The reaction 
xnixture (IS mL final volume) consisted of 1.0 mM 
each of DL-[^^C]pT0lme (M. Kainosho, uiq)ublishcd) 
and L-[^^N]threonme (Isotec). 1.0 mM each of the 
other 18 amino acids, 6.7 (ig/mL of the pK7-Ras 
plasmid with the T7 promoter and the structural gene 
for die Ras protein (Kigawa et al., 1995). 55 mM 
Hepes-KOH (pH 7^, 1.7 mM DTT, 1.2 mM ATP, 
0.8 mM each of CTP, GTP, and UTP, 80 mM creatine 
phosphate (CP) (Boehringer-Mannheim), 250 jtg/mL 
creatine kinase (CK) (Boehringer-Mannheim), 4.0% 
polyethylene glycol (PEG) 8000 (Sigma), 0.64 mM 
3',5'<ycUc AMP, 68 m L(-)-5-formyl-5,6,7.8- 
teirahydrofolic acid (folinic acid), 175 \LgfmL £ coli 
total tRNA (Boehringer-Mannheim), 210 mM potas- 
sium glutamate, 27.5 mM ammonium aceute, and 
10.7 mM magnesium acetate. These components were 
mixed and preincubated at 37 for 5min, after 
which the S30 extract (300 \iUmL) and T7 RNA 
polymerase (93 \ig/mL) were added. The reaction 
was incubated it 37 *C for I h, and was then cen* 
trifuged for 5 min ai 12000 x g. The supernatant 
was desalted on a PD^IO gel-filtradon column (Pbar* 
macia). The Ras protein was purified using an FPLC 
Mono Q anion exchange colunm (Phannacia) fol- 
lowed by an FPLC Superdex 75 gel-filtrarion col- 
unm (I%annacia). The yield of die Ras protein was 
detennined &pm the specific absoibance: A28O = 
0.56 cm~^ mg"^ mL CYamasaki et aL, 1992). The 
large scale reaction for ["ClAsp/[^^N]Ser-selecuve 
labeling was performed in the same manner, ex- 
cept L-[^^C] aspartate (M. Kainosho, unpublished) 
and L-[^^N]serine (Isotec) were- used in place of DL- 
[^^C'lproline and L-['^N]threonlne. 



PrepararionofiheRqfRBD 

For the Raf RED preparation, we used the £L coli 
strain BL21(DE3), which produces the rat Raf- 
1 RBD (aniino acids ,51^^131) containing the ex« 
era vector derived residues MaSMTGGQQMGRGS 
and KLAAALEHHHHHH as its N- and C-terminal 
residues, respectively. The Raf RBD was purified 
from the lysed cell extract by Ni-lOA agarose Ni- 
chelate afi&nity chromatography (Qiagen) followed by 
CM TOYOPEARL cation exchange chromatography 
(Tosoh). The yield of the Raf RBD was detOTiined by 
measuring the tryptophan absorbance at 280 nm. 



The template DNA for in vitro preparation of the 
suppressor 

We constructed a fragment, with die T7 promoter, 
the gene for the amber suppressor IRNA'^^'cua. and 
a B jM cleavage site, from the. cloned £. coli wild- 
type iRNA'^^2 g&tic by polymerase diain reaction with 
diree synthetic primers. The fragment was ligaied into 
pUCll9 (pTRYl)> and its niicleotide sequence was 
confirmed with an ALFied DNA sequencer (Pharma- 
cia). The pTRYl plasmid DNA (1 mg) was digested in 
a reaction mixture (5 mL) with Bsi^ (1500 units) at 
55 ^C for 17 h, and was phenol-chloroform extracted, 
ohanol-precipitated, and dried. 

In vitro preparanon of the suppressor tRNAJ^^cUA 
The iranscripuon was canied out in a reaction mix- 
ture (10 mL) containing 40 mm Hepes-KOH (pH S.l), 
46 mM MgCl2, 5 mM DTT, 40 ^LgM. BSA, 2 mM 
spermine, 0.01% Triton-XlOO, 6.25 mM each of ATP, 
CTP. GTP, and UTP. 20 mM 5'-GMP. 0.1 mg^ 
BsiNI-digested template DNA. 20 units/mL of RNase 
inhibitor (Takara), and 60 p-g/mL T7 RNA poly- 
merase. Aft^ an incubation at 40 ®C for 5 h, 1 mL of 
500 mM BDTA (pH 8.0) was added, and the reaction 
mixmi^ was desalted on a PD-lO column. The tran- 
scribed tRNA was then purified by DEAE-Toyopearl 
column chromatography (Tosoh) foUowed ^y elec- 
trophoresis on a 20% polyactylamide gel containing 
TMurea. 

Preparation of the E. coli tyrosyl-tRNA synthetase / 
(TyrRS) 

The £ coli tyrS gene was cloned firom the £. coli strain 
HBlOl chromosome by polymerase chain reaction, 
and was ligated into theT7 promoter expression vector 
pT7-7 (pT7-7-iyrS). Tlie K coli Sffain JM109 (DE3) 
(Promega) haibQringpT7-7-lyrS was cultured at 37 *C 
in M9ZB medium (Smdier et aL, 1990) containing 
50 M-g/mL ampicillin. The TyrRS protein was purified 
from die lysed cell extract by DEAE-Sephacel col- 
umn chromatography (Phannacia) followed by FPLC 
Phenyl Superose column chroniatography (Pharma- 
cia). The yield of the purified TyrRS protdn was 
determined by a protein assay (Bio-Rad). 

Small-scale preparation cfthe Tyr^tBNA^^'' qv a 
The aminoacylauon of the £ coli suppressor 
tRNA^OJAi piepared as described above, was car- 
ried out in a reaction mixture (60 m-L) containing 
50 mM Hepes-KOH (pH 7.8), 10 mM Mga2, 4 mM 
AIP. 133 M-M tRNA'^cUA. 280 m-M L-[^*qtyrosine 
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(2.23 GBq/mmol, Moravek), and 2 of die purified 
TyrRS. After an incubation at 37 *C for 30 min, flie 
reaction was stopped by die addition of 12 M^L of a 
solution (pH 4 J) containing 12.4 M acetic add and 
3 M ammonium acetate. The ('^QTyr-tRNA'^cUA, 
logedier wiih die unchaiged tKSpJ^cxjAf was phenol- 
chlorofonn extraaed, and then desalted on a BioSpin- 
6 gel-filtration column (Bio-Rad). The sample was 
divided into seven tubes, ethanol-predpitaied, and 
dried Two tubes were immediately used for the two 
types of small-scale cell-free protein syndiesis exper- 
imems (with unlabeled or ^^C-labeled L-Ieudne, as 
described below). Anodier tube was used to deteimine 
the amount of [^^Cr|Tyr-iRNA'''5^^cuA by liquid sdn- 
tilladon counting of the tridiloroacedc add insoluble 
material. 

The template DNA for site-directed incorporation 
The pK7-Ra5 plasmid was mutated at the Tyr^^ codon 
CIAQ of ihe Ras protein to an amber codon (TAG) 
(pK7-RasY32am). The nucleodde sequence was con- 
firmed with an ALFred DNA sequencer. 

Small'scale cell-free reaction for die^directed 
incorporation 

The pellet of die [^^CJiyr-tRNA'^y^cuA sample in a 
mbe (described above) was dissolved to a final con- 
centration of 40 in a soludon (30 \iL) containing 
55 mM Hcpes-KOH (pH 7-5), 2 mM DTT, 12 mM 
AlP, 0.8 mM each of CTP, GTP, and UTP, 80 mM 
CP. 250 ng/mL CK, 4.0% PEG 8000, 0.64 mM 
3\5'-cyclic AMP, 68 \lU folinic acid, 175 p.g/mL 
£ coli total tRNA, 210 mM potassium glutamate, 
27^ mM ammonium acetate, lOJmMiIiagnesiumac'- ^ 
etate, 0.64 mM unlabeled L-leudne, 1 mM L-tyrosine, 
0.7 mM each of the other 18 amino adds, 6.7 |ig/mL 
of die pK7-Rasy32am plasmid, 133 jig^ T7 RNA 
polymerase, and 7.2 yih S30 extracL 

Tliis leacdon mixmre was incubated at 37 **€ for 
1 h. At 5. 10,20, '30,^45, and 60 min, 3.7 jiL^aliquots of 
die leaiction mixture, containing bodi the synthesized 
protdn labeled wifli L-[^^Cltyrosine and the remaining 
I^[^^C]Tyr-cRNA, were removed, mixed with-60 m^L 
of O.IN NaOH, and kept at 25 **C for 30 min to 
deacylaie die ladioactive aminoacyUtRNA, The in- 
corporanon of L-[^^C]iyrosine into the Ras protdn 
was determined by liquid scindlladon counting of the 
tricbloroacenc acid insoluble material, The amount 
of the [^^C]Tyr-containing Ras protein was estimated 
based on die assumption diat the L-[^^(3tyrosine was 
incorporated only into die posidon of lyr^^. The 
amount of residual [^^C]Tyr-tRNA''^cuA was esti- 



mated &om die difference between the L-[^^C]tyrosine 
incorporations widi and widiout die deacylation treat- 
ment 

The incorpoiadon of L-[^^C]leucine into die Ras 
protein was determined by subtraction of die L- 
[^^C]tyrosine incoiporation from the total incorpo- 
ration of L-[^^C]'tyiosine and L-L^'^CJleucine by the 
same reaction, except for die use of 0.64 mM 
[^^Clleudne (193 MBq/mmoI, Amersham) in place of 
the unlabeled L-leudne in die reaction mixture. The 
amount of die [^^CjLeu-containing Ras protein was 
estimated on the jassumption that the L-[^^C]leucine 
was incorporated ^to all Leu residues. 
~ In addition, die Ras protein was synthesized by die 
cell-free coupled transcripdon-transladon system, in 
die same reaction mixtuie as that of the siie^diiected 
L-[^^C]iyrosine incorporation, except for the use of 
die pK7-Ras plasmid, which lacks the amber mutation, 
in place of pK7-RasY32am, 0.64 mM L-[^^qieucine 
(193 MBq/mmol) in place of die unlabeled L-leucine, 
and no [^^qTyr-tiRNA*^ciiA. 

These cell-free reaction products were also ana- 
lyzed by SDS-EaGE widi iricine buffer (Schagger and 
von Jagow, 1987), followed by autoradiography widi 
a Bio Image Analyzer BAS-2000 system (Fuji Film), 

Large-scale preparation ofl^^NTfyr-tRNAj^'cUA • 
For die syndiesis of die [^^N]Tyr-tRNA"^aJA" 2.5 mL 
of 50 n^ Hepes-KOH buffer (pH 7.8), containing 
10 mM MgQi, 6 mM An> 430 \lM IRNA'^cua. 
and 0.6 \iM TyrRS, was prepared and' saturated widi 
l^['^tyrosine (Icon) (2 mg of L-[^^N]^sine pow- 
der was added to the buffer, but was not completely 
dissolved). The mixnire was incubated at 37 °C for 
30 min, and die reaction was stopped by die addition 
of 50 \iL of a solution (pH 4^) containing 12.4 M 
acetic acid and 3 M ammonium acetate. The mixture 
was phenol-chloroformextracted, desalted on a FD-10 
column, ethanol-predpitaied, and dried. The amount 
of syndiesizcd Tyr-KNA'^cUA was determined by 
die same aminoaqylapon reaction in a 50-(lL volume, 
except for die usq of 0.92 mM L-[^*(3tyiosine (1.78 
GBq/mmol) in place of L-[^^N]^sine, followed by 
liquid scintillation counting of the trichloroacetic acid 
insoluble material.. 

Preparation off^^Nftyr^^-Ras protein using the 
cell-free protein synthesis system 
The [^^N]Tyr-tRNA^cuA pellet, prepared as de- 
scribed above, was dissolved to a final concentration 
of 36 \iM in a sdlution (30 mL) containing 55 mM 
Hepes-KOH (pH 75), 2 mM DTT, 1.2 mM ATP. 
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0.8 mM each of CTP, OTP, and UTP, 80 mM CP, 
250 jigtoL CK» 4.0% PEG 8000, 0.64 mM 3'^'- 
cydic AMP, 68 |iM folinic acid, 175 jtg/mL E. coli 
total tRNA, 210 xnM potassium glutaiDaie, 27 J mM 
animomum acetate, 10.7 mM magnesium acetate, 
10,8 mg L-tyrosine, 0.7 mM each of ihc other 19 
amino acids, 6-7 \iglmL of the pK7-RasY32am plas- 
mid, 133 (ig/mLT7RNA polymerase, and 7^ mLS30 
extract. Ilus reaction mixture was incubated at 37 °C 
for 30 TTi^n, and was desalted on a PD-10 column. 
The labeled Ras protein was purified to homogeneity 
with Resource Q (Pharmacia) and S upodex 75 FPLC 
columns. Samples from each purification step were 
analyzed by SDS-PAGE. The yield of the [^^NJiyt^^- 
Ras protein was estimated to be 2.2 mg by the protein 
assay. The unlabeled Ras protein was prepared as 
described previously (Ha ei al., 1989). 

Mass spearomexry 

The purified Ras proteins were desalted by ultrafiltra- 
tion using a Centricon-10 unit (Amicon; Mr cutoff = 
1 X IC)^) as described (Akashi et al., 1996). The ma- 
trix assisted laser desorption ionization-time of flight 
(MALDI-TOF) mass spectra were then measured and 
analyzed on a Bruker REFUK mass spectrometer 
with a delayed extraction system. The Ras protein 
samples were analyzed in the linear mode with 3^- 
dimethoxy-4-hydn)xy-dnnamic acid as die matrix. 

In order to determine the level of ^^N enrichment 
of the [^^C]Asp7[^^N]Ser selectively labeled Ras pro- 
tein, the protein was digested with trypsin- One of 
the digested fcagmenis consisted of residues 136-148, 
containing two Ser residues (Ser^^* and Ser^^^) and 
no Asp residue. This fiagmeni was aimlyzed in the 
reflector mode with 2«5-dihydroxy-benzoic acid as the 
matrix. The simulation of the isotopic distribution pat- 
tern for Ra5(i36-148) was cairied out with tlie Bn^ ^ 
Xmass program on a SUN SS5 workstation. The pat- 
terns for the various ^^N contents were calculated 
as linear combinations of the simulated partes for 
the unlabeled Ras(I36-148) and the labeled Ras(136- 
148) wifli two [^^N]Scr residues. To estimate the ^^N 
content of Tyr^^ of flie [^^N]Tyr^.Ras protein, die 
protein was digested with Achromobacter fydcus pro- 
tease I (lysylendopeptidase). The digest containing die 
Ras fragment consisting of residues 17-42, Ras(17- 
42), was analyzed as described above. The simulation 
of the isotopic distribudon pattern for Ras(17-42) was 
carried out as described above. 



NMR samples 

For NMR measurements of selectively ^^C/*^N- 
labeled Ras proteins, 6 mg of each of the proteins, 
prepared as described above. «^ used To prepare the 
guanosine 5'-0-(P,v-imidotriphosphate)-bound form 
of the Ras protein (Ras«GI^PNP; GMPPNP = 
a slowly hydrolyzable analog of GTP), die bound 
GDP was exchanged by incubating the Ras protein 
(OJ mM) diree times with GMPPNP (Boehringer- 
Mannheim, final concentration!! mM) in the presence 
of 10 mM EDTA and 10 uniti/mL apyrase (Sigma) 
for 10 min at 37 'C. The protein samples were 
dien concentrated and mixed v^tb 6 mg of Raf RED 
(Ras#GMPPNP: Raf RBD 1:2) in NMR buffer 
(90% H2O and 10% D2O coniBining-20 mM sodium 
phosphate buflfer (pH 6.5), 5 mM DTT, 5 mM MgOa, 
100 mM NaCl and 0.01% NaNj) using a Centricon-3 
unit (Amicon; Mr cutoff = 3 x 10^). 

For an NMR measurement of [^^N]Tyr^^-Ras pro- 
lein, 2.2 mg of the protein was used. The solution of 
the Ras protein (0.2 mL) was mixed witii 5 mL of 
the NMR buflfer described above. The solution was 
concentrated to about 0.2 mL by ultrafiltration using 
a Cenlricon-10 unit (Amicon; Mr cutoff = 1 x 10*). 
This cycle was repeated five^times. 

/ 

NMR measurements J 

Spectra of selectively ^^CV^^N-labeled Ras proteins 
w^ recorded on a Bruk^ AMX-^00 spectrome- 
ler. and the spectrum of [^^NJTyr^-Ras protein was 
recorded on a Bruker DMX-500 spectrometer a; a 
probe temperanire of 30 **C. The NMR daia were 
processed and analyzed with die Felix (BiosyxDrMSI), 
the AzaiffCW. Boucher; unpublished), die NMRPipe 
(Delaglio el aL, 1995), and die NMRVIEW (Johnson 
and Blevins, 1994) programs on an Indigo^ worksta- 
tion (SGI). Tlie 2D ^H-^^N HSQC (Bod^ihausen and 
Ruben, 1980) spectra were acquired with 64 (rj) x 
1024 (ri) complex points using the water flyback 
scheme (Grzesiek and Bax. 1993). The 2D ^H-^^N 
HNCO (Ikura et al., 1990) spectra were acquired 
with 32 (ri) x 512 (zi) complex points. For each 
FID, 256 transients for HS<2C, and 1024 transients 
for HNCO were accumulated. By zero-filling followed 
by Fourier transformation, spectra of 1024 (^H) x 
128 (^f N) points for selectively ^^C/^^N-labeled Ras 
proteins and of 1024 (*H) x '512 (^^N) poi^ts for 
[^^•N]Tyr^^-Ras protein, respectively, were obtained. 
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Figure L Iht Uotopic distribution paoems for (he Ra&( 136-1 48) 
fr^gmenL (A) The observed ^ern of the MaLDI-TOF 
mass ^peccmm for the Ras<136-148) fragirtCQt of the 
[^C]Aqi/C^NlSer-»decdveIy labeled Ras j«otciji. Simulaied 
panems for the Ras(l36-148) BagmcDt with various ^^N comenis: 
0% (B), 25% (Q. 50% P). 759fc (E). and 100% (F). 

Results and Discusaon 

Preparation of the selectively ^^C/^^N-labeled Ras 
protein 

For the [^^C]PrQ/[^^N]Thr-selecuve labeling of the 
Ras protein, a IS-mL reaction of the 'final' system 
was used with 2 mg each of DL-[^^C]proline and t- 
[*^N] threonine. Tlie labeled Ras protein was purified 
to homogeaeiQr from the crude cdl-fcee reaction mix- 
tuxe by successive chromaiogiapby on Mono Q and 
Superdex 75 columns. Finally, we obtained approx- 
imately 6 mg of the Ras protein, as determined by 
absorb^nce at 2S0 nzn (data not shown). 

Similarly, for die ['^C]Asp/[^^N]Ser-selective la- 
beling,; a IS-mL* reaction with 2 mg each of L- 
[^^GTflspartate and L-[^^N]serine was used. We ob- 
tained about 6 mg of the Ras protein (data not shown). 

Mass spectrum of the [^^CjAsp/f^^NJSer-selectively 
labeled Ras protein 

To determine the level of ^^N enrichment of the 
[^^C] Asp/[^^N]Ser-selec£ively labeled Ras protein, an 
aliquot of (he sample was digested by oypsin, and was 
subjected to MAlDI-TOF mass spectrometry. The 
isotopic distribution patterns of the peak for the frag- 
ment including two residues and no Asp residue, 



Ras(136-I48), were; analyzed (Figure 1 A). As the iso- 
topic disoibution pa[ttem depends on the ^^N content, 
it is possible to estimate the acmal ^^N content from 
the observed pan^'by comparison with'tbe<simulated 
paneras from various ^^N conieats (Figures IB-F). 
The observed pattern firom the labeled Ras(136-14S) 
was most similar to die simulated pattern for the 75%- 
labeled Ras(136-l48) (Figure IE). Thus, the level 
of ^^N enrichment was estimated' to be 75%. This 
enrichment level is high^ than the case of the dual 
labeling by an in vivo system, in which die level of ^^N 
enrichment was less than 50% (Kainosho and Tsuji, 

1982). ; 

[^^CJPro/[^^N]Thr^electivefy labeled Ras protein 
To measure die ^H^^^N HSQC and ^H-^^N HNCO 
spectra (Figures 2A*and B, respectively), 6 mg each 
of die labeled Ras«GMPPN? and the Raf RBD were 
used (the molar lado was about 1 to 2). In the case 
of the Ras protein produced with L-[^^N]thxeonine by 
the £. coli expression system, die Gly and Ser residues 
were labeled in addition to die Thr residues, because 
the direonine aldolase cataly:^ the durecnine-gtycine 
conversion and theiserine hydroxymethyltransferase 
catalyzes die glycine-serine conversion (Mcintosh 
andDahlquisL 1990; Yamasaldetal., 1992). The num- 
b^ of observed crqss pealcs in the^HSQC spectrum 
(Rgure 2A) is equal co that of the Thr residues in 
the Ras protein, in^cadng diat these metabolic con- 
versions do not occar in the cell-free system. As the 
Pro-Thr sequence lis unique in die Ras ^protein, the 
Thr^^ cross peak was unambiguously identified by the 
HNCO experiment (ipigurc 2B). | 

[^^C]Asp/[^^NJSer'Selectivefy labeled Ras protein 
Tlie *H-"N HSCJC ,and ^H-'^N HNCO spectra (Fig- 
ures 3A and B, respectively) were recorded using 6 mg 
each of die labeled Ras^GMPFNP and die Raf RBD 
(die molar ratio was about 1 to 2). Although the ^^N 
enridiment level of die /sample is sufficiendy high 
(about 75%), die signal-io-noise ratio of die HSQC 
spectmm was not as high as in ibe spectrum of die Ras 
protein by itself at the same Ras concentration (about 
1 mM) (Muto et aL, 1993; Ito et al., 1997). This is 
because of die long correlation time (approximately 
21 ns; Teiada et al., ippublished) of die Ras«Raf RBD 
complex (die total molecular mass is about 30 kDa). 
In the HS(2C speccpim, seven peaks are observed, 
while there are eight serine residues in the Ras pro- 
tein (Figure 3A). Since isotopic transfer from serine 
to any odier amino add could not occur in the cell- 
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Ft^ure 2. The HSQC spcccmm (A) and the HNCO spectrum (B) of r"c']Pro4^*NXni«elea^ labeled Ras^GMPPNP bound io Raf RBD. 



free syscem (Kigawa et al., 1995), we concluded diat 
the seven peaks were derived from tbe Ser residues, 
and two of them weie overlapped. There are two Asp- 
Ser sequences in ihe Ras prorein: Asp^*-Ser^ and 
Asp^^-Ser^^, and two cross peaks were observed in 
die HNCO spectrum (Figure 3B). The peak at ^HN 
= 7-51 ppm and = 110.64 ppm was assigned to 
Ser**^, by comparison with the assignment of Sw*^ 
in both the GDP-bound Ras protein (Ras#GDP, Maro 
et aU 1993; Ito et al., 1997) and Ras»GMPPNP Qxo 
et aL» 1997). This is also supported by tbe considera- 
tion that Ser^^is not located on the binding interface 
with Raf RBD, from tbe crystal structure of the com- 
plex of GMPPNP-bound RaplA (RaplA^GMPPNP), 
which is a homologue of Ras, and the Raf RBD (Na$^ 
sar e: aL, 1995, 1996). The peak at ^HN ^ 9.03 ppm 
and '^N = II 8.04 ppm was assigned to Sor^^, 

Effector region of the Ras proiein complexed mth the 
Raf RBD 

The backbone amide res'onanoes of die residues in 
the effector region (including Thr^^ and-Ser*^)and 
some other regions of Ras*GMPPNP were extremely 
broadened (Ito et al.. 1997). This is because these re- 
gions slowly interconvert between two or more stable 
conformers (^regional polysterisra% Ito et al., 1997). 
On the other hand, three cross peaks, Asp^^, Asp^^, 
arid Asp^^, which were too broad to be detected 
for Ras»GMPPNP by itself (Ito et al., 1997), were 
clearly observed for Ras#GMPPNP in the compile 
with Raf RBD. Moreover, the Thr^^ and Ser^* cross 



peaks could also be observed, as demonstrated in this 
smdy. Thus, tbe conformations of the effiector region 
of RasvGMPPNP are mostly fixed upon binding with 
RafRBD. ; * • 

The diemical shift difEerences of Tbr^ and 
Sei^^between Ras^GMITNP complexed wi^ the Raf 
RBD and Ras#GDP (Muto ef aL, 1993; Ilo ec al^ 
1997) wCTe large ata^^ A^HN « 0.97 ppm and A^^N 
= 16.15 ppm; Ser^, A^HN = 0.67 ppm and A^^N 
« 2.40 ppm), while that of iSer^*^ was negligibly 
small (A^HN = 0.07 ppm arid A% = 0.08 ppm). 
This mdifay^ that the em^ironbents of the Thi^^ and 
Ser^^ amide groups dramatically change upon bind- 
ing widi Raf RBD. In the tertiary structure froni Ae 
X-ray crystallography study, the Thr^^ amide group 
of Ras^JjMPPNP fbnns a hydtpgea bond wifli the y- 
phosphate group of GMPPNP ([Pai et aL, 1989X while 
that of Ras^GDP does not interact with any atoms 
(Nfilbum et aL, 1990). This, is also ihe case with 
RaplA»GMPENP complexed' with Raf RBD (Nas- 
. sar et aL, 1995, 1996). Our result suggests Oiat the 
Thr^^ of Ras«GMPPNP complexed wifli Raf RBD 
also forms a hydrogen bond with the y-phosphate 
group of the GMPPNP. The orientation of Sex^^ of 
the Ras protein is abnost the same between the GDP- 
bound and GMPPNP-bound fprms, according to die 
X-iay smdy (Pai ei al., 1989; NDIbum ei aL, 1990). In 
the RaplA^GMPPNP structure in the complex with 
Raf RBD, the amide group of Ser^^ of RaplA forms 
a hydrogen bond with the caibonyl group of Arg^ 
of Raf RBD (Nassar et aL» 1995, 1996). Our result 
suggests that, in Ras^GMPFNP complexed with Raf 
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rt^ureS.Tht HSQC specmim (A) and ibe HNCO specuum (B) of [^^C]Asp/['^N]SefKselecuvely labeled Ras*GMPFNP bound to Raf RBD. 



RBD, Ser^^ also forms a hydrogen boDd with some 
residue (presumably Aig^') of Raf RBD. 

Scheme for the site'directed ^^N-labeling ofTyr^^ qf 
the Ras protein 

Figure 4 summarizes die procedoies for die site- 
direcied siable-isoiope amino acid labeling of che 
Ras protein. The codon for Tyx^ of the Ras pro- 
tein was replaced by an amber tennination codon 
(TAG) on the plasmid template (pK7-RasY32am). An 
£. coli tyrosine suppressor, iRNA^^cuA* was pre- 
pared by in vitro iranscripdon widi T7 RNA poly- 
merase (25 mg purified tRNA^cuA per 10-mL tran- 
scripiion reaction mixture)- The iRNA'^cua was 
precharged widi L-[^^N]tyiosine by die£ coUTy^, 
The cell-free coupled transciipdon-traAslation reac* 
don was carried out widi [^^NJTyr-tRNA^y'cUA 
pK7-Ra5Y32am zs die ten^late DMA; die amber 
codon (UAG) in die mRNA was recognized by die 
[^^Tyr-iRNA^CUA,and diea L-t^^fjjq^sjjjg ^ 
site-specifically- incorporated into Tyi^^ of die Ras 
protein. The amber codon was conq>edtivdy recog- 
nized by release factor 1 (RFl), and the abordve 
product, the Ras fragment consisting of residues 1« 
31, Ra5(l-31), was also synthesized. The details of 
the cransladon experiments will be described below. 

Time courses ofcell-firee protein synthesis 
In order to opdmize die xeacdon time for the sice- 
diiected isotope labeling (SDIL), die translation ki- 
netics were ^lalyzed on a small scale by the use 
of L-[^^C]iyrosine and L-[^*C]leucine. instead of L- 



[^^N]tyrosine and unlabeled Weucine, respecdvely, 
as in die following. More dian S0% of die suppressor 
tRNA^OJA was aminoacylated widi L-[^^C]tyrosine, 
according to die small-scab preparation protocol (see 
Materials and Methods). The [^^Cjiyr-coniaining pro- 
tein was efficiendy syndiesized widi die SDIL sys- 
tem; the accumulated amount increased, and finally 
reached a plateau of about 100 ^.g/mL at about 45 min 
(Hgure 5, crosses). This amoiini is approximately 
one-third of diai syndiesized in 45 min with the con- 
ventional synthesis system with no amber mutation 
(Figure 5, filled 'circles). Although d^e syndiesis of 
the ['^qTyr-comaining protein stopped at about 45 
min, die syndiesis of the [^^C]Leu-contaiiung pro- 
tein widi the SDJL system sdll continued (Figure 5, 
open circles). At 45 min, die [^*qTyr-tRNA'^cUA 
appears to be exhausted (Figure 5, filled squares), 
and concomitandy die syndiesis of die [^^C]Tyr- 
containing protein stopped. Th^efore* die content of 
die ['^CJiyr-containing protein was relatively high 
(about 70%) until 30 i^in, but then decreased because 
the unlabeled protein was syndiesized by recbazging of 
die deacylated suppressor tRNA'^^'cuA with unlabeled 
L-tyrosine by the endogenous TyrRS. 

SDS'PAGE analysis of the cell-free synthesized 
proteins 

The proteins synthesized widi I^'^ClTyr-iRNA^'cuA 
and/or L-[^^C]leucine by die small-scale reaction 
protocol were analyzed by SDS-PACX with tricine 
buffer, followed by autoradiography (Figure 6). When 
die cell-free protein syndiesis with pK7-RasY32am 



303 



Texnplue DNA for TempUte DNA 

' JtAG suppressor iRNaJJ^ 



r ^ 

lpK7-RasY32am| 
|t7I 



in vtr/o nanscnp&oQ 
by T7 RNA pdymcFftse 



RNApdyxneiase 
eodoni CDdoA32 
— UA C U AG" 
mRNA 



codon 172 
~UAA«» 



^^^^ ^^^^ 



Tyr-iRNASgy, 



Release factor 



r £ coJi cell-free protein synthesis . 




Aminoacyladoa 



t^^NIiyr^lNASU 



( coupled ccai^paoiKn^huioo) 




[^N]TV^'2-RaipfOlfim 
Figure 4. Site-diiected Ktable^sotope l&bdiDf of the Ras pfOtc^ 



lacked [^*C]Tyr-tRNA'^cUA. only the abonive prod- 
uct, Ras(l-31), was syntbesiwd (Figure 6, lane 3). In 
coQtrasL the Ras protein was primarily produced in 
the presence of [^^]Tyr-tRNA'^cuA (Figure 6, lane 
4). The mobility of this Ras protein was confinned to 
be equal to that of the authentic Ras protein, which 
was produced with the conventional synthesis system 
containing the original pK7-Ras plasmid without the^ 
amber mutation (Figure 6, lane 2). The suppression 
efficiency was estimated to be. higher than 90% by 
analysis of die autoradiogram; [^^CJTyi-tRNA^'cUA 
efficiendy suppressed the amber mutadon. It was 
shown that the Ras protein, as synthesized through 
amber suppression, definite!/ contained the [^^C]Tyr 
residue that had been precharged to the suppressor 
tRNA^CUA (Figure 6. lane 5). 

Preparation of the [^^NjTyr^^'Ras protein 

Hie tRNA'^ojA aminoacylated.wiih L-[^^N]iyrosine 
was prepared according to the large-scale prq>aration 
protocol as described in Materials and Methods. As 
the tRNA^ctlA concentration was higher dian that 
in the small-scale reaction protocol, we analyzed the 
aminoacylaiion efficiency firoro die L-['^C]tyiosine in- 



corporation, * and confirmed that about 50%^ pf the 
suppressor iRNA^cua ^as aminoacylaied (data not 
shown). The protein syndiesis reaction was caiiied out 
with [^^NJTyr-iRNA'^cuA- As much as 2.2 mg of the 
Ras protein could be purified fiom a 30-n)L leacdon: 
The small-scale analysis had predicted the produc- 
tion of about 3 mg of crude Ras protein per 30 mL 
reaction, ^nsideting that die 2.2 mg were obtained 
after the protein purification, the present large-scale 
pr^aration of the site-specifically labeled Ras pro- 
tein was quite successful The SDS-PAGE analysis 
revealed diat the purified Ras protdn formed a sin- 
gle band, and its molecular mass was equal to that 
of the authentic Ras protein (Figure 7). This was also 
confirmed by mass spectrometry analysis (data not 
shown). This shows thai the Ras protein was clearly 
separated from die abonive product, Ras(l-3 1). if any» 
by the chromatographic purification procedure. 

Mass specirum of the [^^N]Tyf^^'Ras protein 
In order lo estimate the ^^N content of die Tyr^^ 
residue, labeled and unlabeled Ras proidn samples 
were digested by Achromobacter lyricus protease I, 
and were subjected to MALDI-TOF mass spcctrom- 
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figure S, Tixne coun^es of cdl-£ree protein synthesis. Bach 
of synthesized protem at a given dxne is normalized lo i^g 
Ras protdn per I mL of the xeaction mixaire Qefi-band ordi- 
n&tc), and is al£o indicsued as \M units (right-hand otdio&ie). 
The amount of ('^QTyr-diNA^cuA « normalized to »iM. The 
[^^QTyr-ooniaimng protein (x) and the tool piOKin (O) (see Ma- 
terials axid Methods) as synthesized ^ die cdl-^ protein synthesis 
system involving the amber-muxsni pIC7-RasY32am pUsmid as 
^ ceo^laie toeeiher with jbe [^*C]Tyr-tRNA^cUA SOIL 
:fyEtem), the protein as synthesized by the convemional synthesis 
system with (he pK7-IUs plasmid wiibouc the amber mutation (•), 
and the amount of (^^C]Tyr-iRNA^cUA >^ SDE. system (■). 




+ + 
+ + 
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pK7-Ras ' • + ' 
pK7-RafiY32am - - H 
[14qrryi*NA2[^ . . ■ 

L-[t'4p]teudne + 19 other ^ ^ ^ 
amfnoadds 
20 amino adds 

Figure 6. T^cin&^D$^AG£ analysis of the cell-£ree leacooa 
products (autoradiogFaphy). The componeacs in the reactions axe 
indicaTfd for eadi sampla 



etry. The isocopic disoibudon patterns of the peak for 
the fragment including Tyr^^, Ras(17-42)» were dien 
analyzed (Figure 8A). The observed pattern from the 
labeled Ras protein was tnosc similar to the simulated 




Figure 7. SDS-PACS analysis of san^les containing the cdl'^free 
synthesized Ras proceia at each purificadon step (CBB soiniog). 
T jinft 1: oude reaction producr, lane 2: Resource Q column cfaro- 
maiography fraction; Uine 3: Supcrdex 75 column chromatography 
ibcdon; 1^ 4: Ae authentic Ras protein that was overexpressed 
in vivo and purified; lane M: molecular mass maikers. 
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Figure 8, T1)e isotopic distribution paaems far the Ras( 17-^2) firag- 
mem induding T^c". (A) The observed pattern of the MALDt-TOF 
mass spectnmi for the Ras(17-42) £tagmetw of the [^^N]Tyr^--Ras 
protein. Simulated panems for die Ras(1 7-^2) fragment with vari- 
ous oontents: 0% (B). 25% (Q. 50% (D), 75% 3»d 100% 
CF). 



partem for the 50%-labeled Ras(17-42) (Figure 8D). 
Thus, the ^^N content was estimated to be about 50%, 
which is slightly lower than thai predicted by the 
small-scale analysis described above (about 70%). We 
are planning to sqlve this problem by lowering the rate 
of non-specific deacylation of [^^K|Tyr-tRNA'^y^cuA, . 
for example, by using overproduced elongation factor 
Tu. Tliese efforts may also reduce the total amotuit of 
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Ft^urt 9. The HSQC spectrum of the I^Njryr^--Ras protein (A) 
and its cross sections (B). Axzou^s indlcaie (be poi^iiioas vtoe we 
have observed, for ihe uiufonnly ^^N-labelcd Ras proieuL cro£ii 
peaks of 7>r residues (Muto et al., 1993: Ito et al., 1997). 



suppressor tRNA required for die site-directed label- 
ing. 

NMR spectrum of the [^^NJTyi^^-Ras pwtein 
The [^^N]Tyr^^-Ras protein (2.2 mg) in the GDP- 
bound form was subjected lo NMR analysis. Tbie 
HSQC spectnim has only one c£oss peak (Rgure 9a). 
The and cliemical shifts of the cross. peak 
were equal to those for the Tyr*^ residue of Ras»GDP 
(Muto et al., 1993: Iio ei al., 1997). No cross peak 
corresponding to Tyr residues other than ly^^ was 
observed in the spectrum (Figure 9B). Thus, lyr^^ 
was labeled site-specifically. Isotope scrambling of L- 
[^^NJiyrosine had been expected. However, evai when 
a number of L-[^^N]tyrosine molecules were released 
from [^^NJTyr-tRNA^>^cuA by hydrolysis prior to in- 
corporation into the protein, they were not recycled 
by the endogenous T/rRS, because of dilution by the 



excess amount of unlabeled L-tyrosine in the cell- 
free protein syndiesis system. Thus, we reduced the 
scrambling of [^^Njtyrosine to a negligible level 

Cell-free protein syrUhesis system for i^MR structure 
analysis 

Hie dual amino add-selecdvc '^C-^^N labeling by 
cell-free protein synthesis is usefiil to decrease the 
number of cross peaks* and ^funhennore for reliable 
resonance assignment. The present soidy showed that 
diese advantages are much Enhanced by the use of 
cell-free protein synthesis as ^compared with die con- 
ventional recombinant tedbniques. The cell-free dual 
labeling is particularly powerful for strucniral investi- 
gation of a large protein or a piotein complex focusing 
on a residue of functional importance. 

In addition, we succeeded in obtaining a mil- 
ligram quantity of a Ras protein Aat was isotopically 
labeled at a specific residue. The site-directed incor- 
poration of a stable-isotope labeled amino acid other 
than tyrosine is possible if the'aminoacylated suppres- 
sor tRNA is available. In this context, a number of 
stqqiressor tRNAs, which are effective in vivo, have 
b^ reported (Kleina et aL, 1990; Normanly et al., 
1990). Some of them show low aminoacylation ac- 
tivities in vitro (Shulman, 1991)- For amino adds for 
which no efficient suppressor tRNA is available, the 
aminoacylaied suppressor iRNA could be obtained by 
the chemical aminoacylation ^method (Ellman et al., 
1991; Robertson et aL, 1991; Lodder et aL, 1997) 
or by aminoacylation with a ciiutant aminoacyl-tiE(NA 
syndietase engineered to aminoacylate die suppressor 
tRNA (Schmitt et aL, 1993; Liu et aL, 1997). Sixe- 
directed stable-isotope labeling drastically simplifies 
the bbservauon and resonance^assignment procedures 
foraspecified amino acidresidueof pazticularinterest. 
and will therefore be usefril, for example, for ana- 
lyzing local sirucmres of large proteins and protein- 
protein interactions. Site-directed labeling win also be 
useful for some solid-state NMR methods, such as ac- 
curate analyses of the distance between two specified 
atoms (Peersen et aL, 1992) and'of the oriwitadonal 
constraint of a defined bond vector in an oriented 
sample (Ketchem et al., 1996). Thus, dual labeling 
and site-direcied labeling by cell-free protein synthesis 
will be useful techniques for analyzing the structures 
of proteins. 
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Condusians 

Dual selective ^^C-^^N labeling of the Ras proiein 
was adiievcd using cell-free proiein synthesis. Tliis 
allowed us to investigate the stnictuie of the effec- 
tor region of die Ras protein complexed with Raf 
R£D. We have also synthesized milligraTn quaniities 
of a site^specifically stable-isotope labeled protein and 
have measured the ^H-^^N HSQC spectrum of the 
protein. This method will be applicable to a wide range 
of solution and solid-state NMR studies. 
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The present invention provides a process for identifying compounds which bind to a specific target molecule. The process includes 
the steps of: a) generating a first two-dimensional '^n/ih NMR correlation spectrum of a '■''N-tabcled target molecule; b) exposing the 
labeled target molecule to one or a mixture of chemical compounds; c) generating a second two-dimensional '^N/'H NMR correlation 
spectrum of the labeled target molecule that has been exposed to one or a mixture of compounds in step (b); and d) comparing said first and 
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indcntifying the presence of one or more compounds that are ligands which have bound to the target molecule. 
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Use of Nuclear Magnetic Resonance to Identify 
Ligands to Target Biomolecules 

Technical Field of the Invention 
The present invention pertains to a method for the screening of compounds 
for biological activity and to the determination of binding dissociation constants 
using two-dimensional ^^N/^H NMR correlation spectral analysis to identify and 
design ligands that bind to a target biomolecule. 

Background of the Invention 

One of the most powerful tools for discovering new drug leads is random 
screening of synthetic chemical and natural product databases to discover 
compounds that bind to a particular target molecule (i.e., the identification of 
ligands of that target). Using this method, ligands may be identified by their ability 
to form a physical association with a target molecule or by their ability to alter a 
function of a target molecule. 

When physical binding is sought, a target molecule is typically exposed to 
one or more compounds suspected of being ligands and assays are performed to 
determine if complexes between the target molecule and one or more of those 
compounds are formed. Such assays, as is well known in the art, test for gross 
changes in the target molecule (e.g., changes in size, charge, mobility) that indicate 
complex formation. 

Where functional changes are measured, assay conditions are established 
that allow for measurement of a biological or chemical event related to the target 
molecule (e.g., enzyme catalyzed reaction, receptor- mediated enzyme activation). 
To identify an alteration, the function of the target molecule is determined before 
and after exposure to the test compounds. 

Existing physical and functional assays have been used successfully to 
identify new drug leads for use in designing therapeutic compounds. There are. 
however, limitations inherent to those assays that compromise their accuracy, 
reliability and efficiency. 

A major shortcoming of existing assays relates to the problem of "false 
positives'*. In a typical fimctional assay, a "false positive" is a compound that 
triggers the assay but which compound is not effective in eliciting the desired 
physiological response. In a typical physical assay, a "false positive" is a 
compound that, for example, attaches itself to the target but in a non-specific 
manner (e.g., non-specific binding). False positives are particularly prevalent and 
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problematic when screening higher concenlrations of putative ligands because 
many compounds have non-specific affects at those concentrations. 

[n a similar fashion, existing assays are plagued by the problem of "false 
negatives", which result when a compound gives a negaUve response in the assay 
5 but which compound is actually a ligand for the target. False negatives typically 
occur in assays that use concentrations of lest compounds that are either too high 
(resulting in toxicity) or too low relative to the binding or dissociation constant of 
the compound to the target 

Another major shortcoming of existing assays is the limited amount of 
10 information provided by the assay itself. While the assay may comecay idenUfy 
compounds that attach to or ehcit a response from the target molecule, those assays 
typically do not provide any information about either specific binding sites on the 
target molecule or structure activity relationships between the compound being 
tested and the target molecule. The inability to provide any such informauon is 
15 particularly problematic where the screening assay is being used to identify leads 
for further study. 

It has recently been suggested that X-ray crystallography can be used to 
identify the binding sites of organic solvents on macromolccules. However, this 
method cannot determine the relative binding affinities at different sites on the 
20 target. It is only applicable to very stable target proteins that do not denature in the 
presence of high concentrations of organic solvents. Moreover, this approach is 
not a screening method for rapidly tesUng many compounds that are chemically 
diverse, but is limited to mapping the binding sites of only a few organic solvents 
due to the long time needed to determine the individual crystal structures. 

Compounds are screened to identify leads that can be used in the design of 
new drugs that alter the function of the target biomolecule. Those new drugs can 
be structural analogs of identified leads or can be conjugates of one or more such 
lead compounds. Because of the problems inherent to existing screening methods, 
those methods are often of little help in designing new drugs. 

There continues to be a need to provide new, rapid, efficient, accurate and 
reliable means of screening compounds to identify and design ligands that 
specifically bind to a particular target 
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Brief Summary of the Invention 
In one aspect, the present invention provides a process of screening 
compounds for biological activity to identify ligands that bind to a specific target 
molecule. That process comprises the steps of: a) generating a first two- 
5 dimensional ^^N/^H NMR correlation spectrum of a ^^N-labeled target molecule; 

b) exposing the labeled target molecule to one or a mixture of chemical compounds; 

c) generating a second two-dimensional ^^N/^H NMR correlation spectrum of the 

labeled tarcet molecule that has been exposed to one or a mixture of compounds in 

15 1 

step (b) ; and d) comparing said first and second two-dimensional N/ H NMR 

10 correlation spectra to determine differences between said first and said second 

spectra, the differences identifying the presence of one or more compounds that are 
ligands which have bound to the target molecule. 

Where the process of the present invention screens more than one 
compound in step (b), that is. a mixture of compounds, and where a difference 

1 5 between the first spectrum generated from the target molecule alone and that 

generated from the target molecule in the presence of the mixture, additional steps 
are performed to identify which specific compound or compounds contained in the 
mixture is binding to the target molecule. Those additional steps comprise the 
steps of e) exposing the l^N-labeled target molecule individually to each compound 

20 of the mixture, 0 generating a two-dimensional ^^N/ ^ H NMR correlation 

spectrum of the labeled target molecule that has been individually exposed to each 
compound; and g) comparing each spectrum generated in step 0 to the first 
spectrum generated from the target molecule alone to determine differences in any 
of those compared spectra, the differences identifying the presence of a compound 

25 that is a licand which has bound to the target molecule. 

15 I / 

Because the chemical shift values of the particular N/ H signals in the 
two-dimensional correlation spectrum correspond to known specific locations of 
atomic groupings in the target molecule (e.g.. the N-H atoms of the amide or 
peptide link of a particular amino acid residue in a polypeptide), the process of the 

30 present invention allows not only for the for identification of which compound(s) 
bind to a particular target molecule,* but also permit the determination of the 
particular binding site of the ligand on the target molecule. 

In a second aspect, the present inventicxi provides a process of determining 
the dissociation constant, Kd, for a given ligand and its target molecule. Thai 

35 process comprises the steps of a) generating a first two-dimensional ^^N/^ H NMR 
correlation spectrum of a ^^N-labelcd target molecule; b) exposing the labeled 
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target molecule to various concentrations of a ligand; c) generating a two- 
dimensional ^^N/^H NMR correlation spectrum at each concentration of ligand in 
step (b); d) comparing each spectrum from step (c) to the first spectrum from step 
(a); and e) calculating the dissociation constant between the target molecule and the 
5 ligand from those differences according to the equation: 

Kn= ([P3o-x)([Llo-x) 

An advantageous aspect of the present invention is the capability of the 
process of the present invention to determine the dissociation constant of one 
10 ligand of the target molecule in the presence of a second molecule already bound to 
the ligand. This is generally not possible with prior ait methods which employ 
"wet chemical" analytical methods of determining binding of a ligand to a target 
mcdecule substrate. 

In this preferred embodiment, the process of determining the dissociation 
15 constant of a ligand can be performed in the presence of a second bound ligand. In 
accordance with this embodiment, the *^N-labeled target molecule is bound to that 
second ligand before exposing that target to the test compounds. 

The ability of the present method to determine not only the existence of 
binding between one hgand and the target molecule, but also the particular site of 
20 binding in the presence of a second bound ligand permits the capability to design a 
drug that comprises two or more linked moieties made up of the ligands. 

This method uses the two-dimensional NMR correlation 

spectroscopic screening process as set forth above to identify a first and 
subsequent ligands that bind to the target molecule. A complex of the target 
25 molecule and two or more ligands is formed and the three-dimensional suucture of 
that complex is determined preferably using NMR spectroscopy or X-ray 
crystallography. That three-dimensional structure is used to determine the spatial 
orientaUon of the ligands relative to each other and to the target molecule. 

Based on the spatial orientation, the ligands are linked together to form the 
30 drug. The selection of an appropriate linking group is made by maintaining the 
spaUal orientation of the ligands to one anolhcr.and to the target molecule based 
upon principles of bond angle and bond length information well known in the 
organic chemical art 
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Thus, the molecular design method comprises identifying a first ligand 
moiety to the target molecule using two-dimensional ^^N/^H NMR correlation 
spectroscopy; identifying subsequent ligand moieties to the target molecule using 
two-dimensional ^^N/^H NMR correlation spectroscopy; forming a complex of 
5 the first and subsequent ligand moieties to the target molecule; determining the 
three dimensional structure of the complex and, thus, the spatial orientation of the 
first and subsequent ligand moieties on the target molecule; and linking the first and 
subsequent ligand moieties to form the drug to maintain the spatial orientation of 
the ligand moieties. 

10 The identification of subsequent ligand moieities can be perfonned in the 

absence or presence of the first ligand (e.g.. the target molecule can be bound to 
the first ligand before being exposed to the test compounds for identification of the 
second ligand). 

In a preferred embodiment, the target molecule used in a screening or 
1 5 design process is a polypeptide. The polypeptide target is preferably produced in 
recombinant form from a host cell transformed with an expression vector that 
contains a polynucleotide that encodes the polypeptide, by culturing tiie 
transformed host cell in a medium that contains an assimilable source of N such 
that the recombinantly produced polypeptide is labeled with N . 

20 

Brief Description of the Drawings 
In the drawings which form a portion of the specification: 
RG. 1 shows a ^^N/^H correlation spectrum of the DNA binding domain of 
uniformly ^^N-labeled human papillomavirus E2. The spectrum (80 
25 complex points, 4 scans/fid) was acquired on a 0.5 mM sample of E2 in 20 

mM phosphate (pH 6.5), 10 mM dilhiothreitol (DTT) and 10% deuterium 
oxide (D2O). 

FIG. 2 shows ^^N/^H correlation spectra of the DNA binding domain of 

uniformly ^^N-labeled human papillomavims E2 before (thin multiple 
30 contours) and after (thick single contours) addition of a final test 

compound Theftnalconcentrationof compound was 1.0 mM. All other 
conditions are as slated in FIG. 1. Selected residues that show significant 
changes upon binding are indicated. 
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FIG. 3 shows ^^N/^H correlation spectra of the DNA binding domain of 

uniformly ^^N-labeled human papillomavirus E2 before (thin multiple 
contours) and after (thick single contours) addition of a second test 
compound. The final concentration of compound was 1.0 mM. All other 
5 conditions are as stated in HG. 1. Selected residues that show significant 

changes upon binding £irc indicated 
FIG. 4 shows ^^N/^H conrelation spectra of the catalytic domain of uniforaily 
^^N-Iabeled stromelysin before (thin multiple contours) and after (thick 
single contours) addition of a test compound. The final concentration of 
10 compound was 1.0 mM. The spectra (80 complex points, 8 scans/fid) 

were acquired on a 0.3 mM sample of SCD in 20 mM TRIS (pH 7.0), 20 
mM CaCl2 and 10% D2O. Selected residues that show significant changes 
upon binding arcindicated. 
FIG. 5 shows ^^N/^H correlation spectra of the Ras-binding domain of uniformly 
15 ^^N-labeled RAF peptide (residues 55-132) before (thin multiple contours) 

and after (thick single contours) addition of a test compound. The final 
concentration of compound was 1.0 mM. The spectra (80 complex points, 
8 scans/fid) were acquired on a 0.3 mM sample of the I^F fragment in 20 
mM phosphate (pH 7,0), 10 mM DTT and 10% D2O. Selected residues 
20 that show significant changes upon binding are indicated. 

FIG. 6 shows ^^N/^H correlation specu-a of uniformly ^^N-labelcdFKBP before 
(thin multiple contours) and after (thick single contours) addition of a test 
compound. The final concentration of compound was 1.0 mM. The 
spectra (80 complex points, 4 scans/fid) was acquired on a 0.3 mM sample 
25 of FKBP in 50 mM phosphate (pH 6.5), 100 mM NaCl and 10% D2O. 

Selected residues that show significant changes upon binding are indicated. 
FIG. 7 shows a first depiction of the NMR-derived smicture of the DNA-binding 
domain of E2. The two monomers of the symmetric dimer are oriented in a 
top-bottom fashion, and the N- and C-termini of each monomer are 
30 indicated (N and C for one monomer. N* and C* for the other). Shown in 

ribbons are the residues which exhibit significant chemical shift changes 
(A6(^H)>0.04 ppm; A6(^^N) >0.1 ppm) upon binding to a first test 
compound. These residues correspond to the DNA-rccognilion helix of 
E2. Selected residues arc numbered for aid m visualization. 
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FIG. 8 shows a second depiction of the NMR-derived structure of the DNA- 
binding domain of E2. The two monomers of the symmetric dimer arc 
oriented in a top-bottom fashion, and the N- and C- termini of each 
monomer are indicated (N and C for one monomer, N* and C* for the 
5 other). Shown in ribbons are the residues which exhibit significant 

chemical shift changes (A6(^H)>0.04 ppm; A6(^^N) >0.1 ppm) upon 
binding to a second test compound. These residues are located primarily in 
the dimer interface region. Selected residues are numbered for aid in 
visualization. 

10 HG. 9 shows a depiction of the NMR-derived structure of the catalytic domain of 
stromelysin. The N- and C-termini are indicated. Shown in ribbons are 
the residues which exhibit significant chemical shift changes (A6(^H)>0.04 
ppm; A6(^^N) >0.1 ppm) upon binding to a test compound. These either 
form part of the ST binding site or are spatially proximal to this site. 

1 5 Selected residues are numbered for aid in visualization. 

FIG. 10 shows a ribbon plot of a ternary complex of first and second ligands 
bound to the catalytic domain of stromelysin. 

Detailed Description of the Invention 
20 The present invention provides a rapid and efficient screening method for 

identifying ligands that bind to therapeutic target molecules. 

Ligands are identified by testing the binding of molecules to a target 
molecule (e.g.. protein, nucleic acid, etc.) by following, with nuclear magnetic 
resonance (NMR) spxjctroscopy, the changes in chemical shifts of the target 
25 molecule upon the addition of the hgand compounds in the database. 

From an analysis of the chemical shift changes of the target molecule as a 
functicMi of ligand concentration, the binding affinities of ligands for biomolecules 
are also determined. 

The location of the binding site for each ligand is determined from an 
30 analysis of the chemical shifts of the biomolecule that change upon the addition of 
the ligand and frcmi nuclear Overhauser effects (NOEs) between the ligand and 
biomolecule. 

Information about the structure/activity relationships between ligands 
identified by such a process can then be used to design new drugs that serve as 
35 ligands to the target molecule. By way of example, where two or more ligands to a 
given target molecule are identified, a complex of those ligands and the target 
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molecule is formed. The spatial orientation of the ligands to each other as well as 
to the target molecule is derived from the three-dimensional structure. That spatial 
orientation defines the distance between the binding sites of the two ligands and the 
orientation of each ligand to those sites. 
5 Using that spatial orientation data, the two or more ligands are then linked 

together to form a new ligand. Linking is accomplished in a manner that maintains 
the spatial orientation of the ligands to one another and to the target molecule. 

There are numerous advantages to the NMR-based discovery process of the 
present invention. First, because a process of the present invenUon identifies 
10 ligands by directly measuring binding to the target molecule, the problem of false 
positives is significantly reduced. Because the present process identifies specific 
binding sites to the target molecule, the problem of false positives resulting from 
the non-specific binding of compounds to the target molecule at high 
concentrations is eliminated. 
15 Second, the problem of false negatives is significanUy reduced because the 

present process can identify compounds that specifically bind to the target molecule 
with a wide range of dissociation constants. The dissociation or binding constant 
for compounds can actually be determined with the present process. 

Other advantages of the present invention result from the variety and 
20 detailed data provided about each ligand from the discovery process. 

Because the location of the bound ligand can be determined from an 
analysis of the chemical shifts of the target molecule that change upon the addiuon 
of the ligand and from nuclear Overhauser effects (NOEs) between the ligand and 
biomolecule, the binding of a second ligand can be measured in the presence of a 
25 first ligand that is already bound to the target. The ability to simultaneously 
identify binding sites of different ligands allows a skilled artisan to 1) define 
negative and positive cooperative binding between ligands and 2) design new 
drugs by linking two or more ligands into a single compound while maintaining a 
proper orientation of the ligands to one another and to their binding sites. 
30 Further, if multiple binding sites exist, the relative affinity of individual 

binding moieties for the different binding sites can be measured from an analysis of 
the chemical shift changes of the target molecule as a function of the added 
concentration of the ligand. By simultaneously screening numerous structural 
analogs of a given compound, detailed sUiicture/aclivily relaUonships about ligands 
35 is provided. 
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In its principal aspect, the present invention provides a process of screening 
compounds to identify ligands that bind to a specific target molecule. That process 
comprises the steps of: a) generating a first two-dimensional ^^N/^H NMR 
correlation spectrum of a ^^N-labeled target molecule; b) exposing the labeled 

5 target molecule to one or more compounds; c) generating a second two- 
dimensional ^^N/^H NMR correlation spectrum of the labeled target molecule that 
has been exposed to the compounds of step (b); and d) comparing the first and 
second spectra to determine whether differences in those two spectra exist, which 
differences indicate the presence of one or more ligands that have bound to the 

10 target molecule. 

Where a process of the present invention screens more than one compound 
in step (b) and where a difference between spectra is observed, additional steps are 
performed to identify which specific compound is binding to the target molecules. 
Those additional steps comprise generating a two-dimensional '^^N/^H NMR 

15 correlation spectrum for each individual compound and comparing each spectnmi 
to the first spectrum to determine whether differences in any of those compared 
spectra exist, which differences indicate the presence of a ligand that has bound to 
the target molecule. 

Any ^^N-Iabeled target molecule can be used in a process of the present 

20 invention. Because of the importance of proteins in medicinal chemistry, a 

preferred target molecule is a polypeptide. The target molecule can be labeled with 
^^N using any means well known in the art. In a preferred embodiment, the target 
molecule is prepared in recombinant form using transformed host cells. In an 
especially preferred embodiment, the target molecule is a polypeptide. Any 

25 polypeptide that gives a high resolultion NMR spectrum and can be partially or 
uniformly labeled with ^^N can be used. The preparation of uniformly ^^N- 
labeled exemplary polypeptide target molecules is set forth hereinafter in the 
Examples. 

A preferred means of preparing adequate quantities of unifonmly ^^N- 
30 labeled polypeptides is to transform a host cell with an expression vector that 

contains a polynucleotide that encodes that polypeptide and culture the transformed 
cell in a culture medium that contains assimilable sources of ^^N, Assimilable 
sources of ^^N are well known in the art. A preferred such source is ^^NH4C1. 
Means for preparing expression vectors that contain polynucleotides 
35 encoding specific polypeptides are well known in the art In a similar maimer, 
means for transforming host cells with those vectors and means for culturing those 
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transformed cells so that the polypeptide is expressed are also well known in the 
art. 

The screening process of the present invention begins with the generation 
or acquisition of a two-dimensional ^^N/^H correlation spectrum of the labeled 
5 target molecule. Means for generating two-dimensional N/ H correlation 
spectra are well known in the art (see. e.g., D. A. Egan et al., BiQchcmisU^, 
32(8): 1920-1927 (1993); Bax, A., Grzesiek, S., Acc. Chcm. Res., 26(4): BI- 
BS (1993)). 

The NMR spectra that are typically recorded in the screening procedure of 

15 1 

1 0 the present invention are two-dimensional N/ H heteronuclear single quantum 
correlation (HSQC) spectra. Because the ^^N/^H signals corresponding to the 
backbone amides of the proteins are usually well-resolved, the chemical shift 
changes for the individual amides are readily monitored. 

In generating such specua, the large water signal is suppressed by spoiling 
1 5 gradients. To facilitate the acquisition of NMR data on a large number of 

compounds (e.g., a database of synthetic or naturally occurring small organic 
compounds), a sample changer is employed. Using the sample changer, a total of 
60 samples can be run unattended. Thus, using the typical acquisition parameters 
(4 scans per free induction decay (fid), 100-120 HSQC spectra can be acquired in a 

20 24 hour period. 

To facilitate processing of the NMR data, computer programs are used to 
transfer and automatically process the multiple two-dimensional NMR data sets, 
including a routine to automatically phase the two-dimensional NMR data. The 
analysis of the data can be facilitated by formatting the data so that ihc individual 
25 HSQC spectra are rapidly viewed and compared to the HSQC spectrum of the 
control sample containing only the vehicle for the added compound (DMSO). but 
no added compound. Detailed descriptions of means of generating such two- 
dimensional ^^N/^H correlation spectra are set forth hereinafter in the Examples. 
A representative two-dimensional ^^N/^H NMR correlation spectrum of an 
30 ^^N-labelcd target molecule (polypeptide) is shown in FIG. I (the DNA-binding 
domain of the E2 protein). 

Following acquisition of the first spectrum, the labeled target molecule is 
exposed to one or more test compounds. Where more than one test compound is 
to be tested simultaneously, it is preferred to use a database of compounds such as 
35 a plurality of small molecules. Such molecules are typically dissolved in 
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perdeuterated dimethylsulfoxide . The compounds in the database can be 
purchased from vendors or created according to desired needs. 

Individual compounds can be selected inter alia on the basis of size 
(molecular weight = 100-300) and molecular diversity. Compounds in the 

5 collection can have different shapes (e.g., flat aromatic rings(s), puckered aliphatic 
rings(s), straight and branched chain aliphalics with single, double, or triple 
bonds) and diverse functional groups (e.g., carboxylic acids, esters, ethers, 
amines, aldehydes, ketones, and various heterocyclic rings) for maximizing the 
possibility of discovering compounds that interact with widely diverse binding 

10 sites. 

The NMR screening process of the present invention utilizes ligand 
concentrations ranging from about 0.1 to about 10.0 mM. At these concentrations, 
compounds which are acidic or basic can significantly change the pH of buffered 
protein solutions. Chemical shifts are sensitive to pH changes as well as direct 

1 5 binding interactions, and "false positive" chemical shift changes, which are not the 
result of ligand binding but of changes in pH, can therefore be observed. It is thus 
necessary to ensure that the pH of the buffered solution does not change upon 
addition of the ligand. One means of controlling pH is set forth below. 

Compounds are stored at 263"K as 1.0 and 0. 1 M stock solutions in 

20 dimethylsulfoxide (DMSO). This is necessary because of the limited solubility of 
the ligands in aqueous solution. It is not possible to directly adjust the pH of the 
DMSO solution. In addition, HCl and NaOH form insoluble salts in DMSO, so 
alternative acids and bases must be used. The following approach has been found 
to result in stable pH. 

25 The 1 .0 M stock solutions in DMSO are diluted 1 : 10 in 50 mM phosphate, 

pH 7.0. The pH of that diluted aliquot solution is measured. If the pH of the 
aliquot is unchanged (i.e., remains at 7.0), a working solution is made by diluting 
the DMSO stock solution 1: 10 to make a 0. 1 M solution and that solution is stored. 

If the pH of the diluted aliquot is less than 7.0, ethanolamine is added to the 
30 1.0 M stock DMSO solution, that stock solution is then diluted 1: 10 with 

phosphate buffer to make another aliquot, and the pH of the aliquot rechecked. 

If the pH of the diluted aliquot is greater than 7,0. acetic acid is added to the 
1.0 M stock DMSO solution, that stock solution is then diluted 1: 10 with 
phosphate buffer to make another aliquot, and the pH of the aliquot rechecked. 
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Ethanolamine and acetic acid are soluble in DMSO. and the proper 
equivalents are added to ensure that upon transfer to aqueous buffer, the pH is 
unchanged. Adjusting the pH is an interactive process, repeated until the desired 
result is obtained. 

5 Note that this procedure is performed on 1: 10 dilutions of 1.0 M stock 

solutions (100 mM ligand) to ensure that no pH changes are observed at the lower 
concentrations used in the experiments (0.1 to 10 mM) or in different/weaker 
buffer systems. 

Following exposure of the ^^N-labeled target molecule to one or more test 
1 0 compounds, a second two-dimensional ^^N/ ^ H NMR correlation spectrum is 
generated. That second spectrum is generated in the same manner as set forth 
above. The first and second spectra are then compared to determine whether there 
are any differences between the two spectra. Differences in the two-dimensional 
^^N/^H NMR correlation spectra that indicate the presence of a ligand correspond 
15 to ^ ^N-labeled sites in the target molecule. Those differences are determined using 
standard procedures well known in the art. 

By way of example, HGs. 2. 3. 4, 5 and 6 show comparisons of 
correlation spectra before and after exposure of various target molecules to various 
test compounds. A detailed description of how these studies were performed can 

20 be found hereinafter in Examples 2 and 3. 

Particular signals in a two-dimensional ^^N/^H correlation spectrum 
correspond to specific nitrogen and proton atoms in the target molecule (e.g., 
particular amides of the amino acid residues in the protein). By way of example, it 
can be seen from FIG. 2 that chemical shifts in a two-dimensional N/ H 
26 correlation of the DNA-binding domain of E2 exposed to a test compound occurred 
at residue positions 15 (115), 21 (Y21), 22 (R22) and 23 (L23). 

It can be seen from HO. 2 that the binding of the ligand involved the 
isoleudne (lie) residue at position 15. the tyrosine (Tyr) residue at position 21. the 
arginine (Arg) residue at position 22 and the leucine (Leu) residue at position 23. 
30 Thus, a process of tiie present invention can also be used to identi fy the specific 
binding site between a ligand and target molecule. 

The region of the protein that is responsible for binding to the individual 
compounds is identified from the particular amide signals that change upon the 
addition of the compounds. These signals are assigned to the individual amide 
35 groups of the protein by standard procedures using a variety of well-established 
heteronuclear multi-dimensional NMR experiments. 
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To discover molecules that bind more tightly to the protein, molecules are 
selected for testing based on the structure/activity relationships from the initial 
screen and/or structural information on the initial leads when bound to the protein. 
By way of example, the initial screening may result in the identification of ligands. 

5 all of which contain an aromatic ring. The second round of screening would then 
use other aromatic molecules as the test compounds. 

As set forth hereinafter in Example 2, an initial screening assay for binding 
to the catalytic domain of stromelysin identified two biaryl compotmds as ligands. 
The second round of screening thus used a series of biaryl derivatives as the test 

10 compounds. 

The second set of test compounds are initially screened at a concentration of 
1 mM. and binding constants arc measured for those that show affinity. Best leads 
that bind to the protein are then compared to the results obtained in a functional 
assay. Those compounds that are suitable leads are chemically modified to 

15 produce analogs with the goal of discovering a new pharmaceutical agent. 

In another aspect, the present invention provides a process for determining 
the dissociation constant between a target molecule and a ligand that binds to that 
target molecule. That process comprises the steps of: a) generating a first two- 
dimensional ^^N/^HNMR correlation spectrum of a ^^N-labelcd target molecule; 

20 b) titrating the labeled target molecule with various concentrations of a ligand; c) 
generating a two-dimensional ^^N/^H NMR correlation spectrum at each 
concentration of ligand from step (b); d) comparing each spectrum from step (c) 
both to the first spectrum from step (a) and to all other spectra from step (c) to 
quantify differences in those spectra as a function of changes in ligand 

25 concentration; and e) calculating the dissociation constant (Kd) between the target 
molecule and the ligand from those differences. 

Because of their importance in medicinal chemistry, a preferred target 
molecule for use in such a process is a polypeptide. In one preferred embodiment, 
a process of determining the dissociation constant of a ligand can be performed in 

30 the presence of a second ligand. In accordance with this embodiment, the '^^N- 
labeled target molecule is bound to that second ligand before exposing that target to 
the test compounds. 

Binding or dissociation constants are measured by following the ^^N/^H 
chemical shifts of the protein as a function of ligand concentration. A known 

35 concentration (P]o) of the target moleule is mixed with a known concentration 
([L]o) of a previously identified ligand and the two-dimensional ^^N/^H 
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correlation spectrum was acquired. From this spectrum, observed chemical shift 
values (6obs) are obtained. The process is repeated for varying concentraUons of 
the ligand to the point of saturation of the target molecule, when possible, in which 
case the limiting chemical shift value for saturation (6sat) is measured. 

5 In those situations where saturation of the target molecule is achieved, the 

dissociation constant for the binding of a particular ligand to the larger molecule is 
calculated using the formula* 

^ _ ([Plo-x) ([L]o-x) 

K.D 

X 

where [PJo is the total molar concentration of target molecule; (L]o is the total molar 
10 concentration of ligand; and x is the molar concentration of the bound species. The 
value of X is determined from the equation: 

x= 



15 where 6free is the chemical shift of the free species; 6obs is the observed chemical 
shift; and A is the difference between the limiting chemical shift value for saturation 
(6sat) and the chemical shift value of the target molecule free of ligand (6free) • 

The dissociation constant is then determined by varying its value until a 
best fit to the observed data is obtained using standard curve-fitting statistical 
20 methods. In the case where 6sat is not directly known, both Kd and 6sai are varied 
and subjected to the same curve-fitting procedure. 

The use of the process of the present invention to detemiine the dissociation 
or binding affinity of various ligands to various target molecules is set forth 
hereinafter in Examples 2 and 3. 
25 Preferred target molecules, means for generating spectra, and means for 

comparing spectra are the same as set forth above. 

The iniUal step in the design process is the identification of two or more 
ligands that bind to the specific target molecule. The identification of such ligands 
is done using two-dimensional ^^N/^H NMR correlation spectroscopy as set forth 
30 above. 

Once two or more ligands are identified as binding to the target molecule at 
different sites, a complex between the target molecule and ligands is formed. 
Where there are two ligands, that complex is a ternary complex. Quaternary and 
other complexes are formed where there are three or more ligands. 
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Complexes are formed by mixing the target molecule simultaneously or 
sequentially with the various ligands under circumstances that allow those ligands 
to bind the target. Means for determining those conditions arc well known in the 
art. 

5 Once that complex is formed, its three-dimensional structure is determined. 

Any means of determining three-dimensional structure can be used. Such methods 
are well known in the art Exemplary and preferred niethods are NMR and X-ray 
crystallography. The use of three-dimensional double- and triple resonance NMR 
to determine the three-dimensional structure of two ligands bound to the catalytic 
10 domain of stromelysin is set forth in detail hereinafter in Example 4. 

An analysis of the three-dimensional structure reveals the spatial orientation 
of the ligands relative to each other as well as to the conformation of the target 
molecule. First, ^he spatial orientation of each ligand to the target molecule allows 
for identification of those portions of the ligand directly involved in binding (i.e., 
15 those portions interacting with the target binding site) and those portions of each 
ligand that project away from the binding site and which portions can be used in 
subsequent linking procedures. 

Second, the spatial orientation data is used to map the positions of each 
ligand relative to each other. In other words, discrete distances between the 
20 spatially oriented ligands can be calculated. 

Third, the spatial orientation data also defines the three-dimensional 
relationships amongst the ligands and the target. Thus, in addition to calculating 
the absolute distances between ligands. the angular orientations of those ligands 
can also be determined. 
25 Knowledge of the spatial orientations of the ligands and target is then used 

to select linkers to link two or more ligands together into a single entity that 
contains all of the ligands. The design of the linkers is based on the distances and 
angular orientation needed to maintain each of the ligand portions of the single 
entity in proper orientation to the target 
30 The three-dimensional conformation of suitable linkers is well known or 

readily ascertainable by one of ordinary skill in the art While it is theoretically 
possible to link two or more ligands together over any range of distance and three- 
dimensional projection, in practice certain limitations of distance and projection are 
preferred. In a preferred embodiment, ligands are separated by a distance of less 
35 than about 15 Angstroms (A), more preferably less than about 10 A and, even 
more preferably less than about 5 A. 
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Once a suitable linker group is identified, the ligands are linked with that 
linker. Means for linking ligands are well known in the art and depend upon the 
chemical structure of the ligand and the linking group itself. Ligands are linked to 
one another using those portions of the ligand not directly involved in binding to 
5 the target molecule. 

A detailed description of the design of a drug that inhibits the proteolytic 
activity of stromelysin, which drug was designed using a process of the present 
invention is set forth hereinafter in Example 4. 

The following Examples illustrate preferred embodiments of the present 
10 invention and are not limiting of the specification and claims in any way. 

Example 1 

Preparatjon Of Uniformly — N-Labeled Targe t Molecules 
A. Stromelvsin 

1 5 Human stromelysin is a 447-amino acid protein believed to be involved in 

proteolytic degradation of cartilage. Cartilage proteolysis is believed to result in 
degradative loss of joint cartilage and the resulting impairment of joint function 
observed in both osteoarthritis and rheumatoid arthritis. The protein possesses a 
series of domains including N-terminal latent and propetide domains, a C-tcrminal 
20 domain homologous with homopexin, and an internal catalytic domain. 

Studies have shown that removal of the N-terminal prosequence of 
approximately eighty amino acids occurs to convert the proenzyme to the 45 kDa 
mature enzyme. Furthermore, studies have shown that the C-terminal homopexin 
homologous domain is not required for proper folding of the catalytic domain or 
25 for interaction with an inhibitor. (See, e.g.. A. I. Marcy, piochcqiistry, 30: 6476- 
6483 (1991). Thus, the 81-256 amino acid residue internal segment of stromelysin 
was selected as the protein fragment for use in identifying compounds which bind 
to and have the potential as acting as inhibitors of stromelysin. 

To employ the method of the present invention, it was necessary to prepare 
30 the 81-256 fragment (SEQ ID NO: 1) of stromelysin in which the peptide backbone 
was isotopically enriched with and ^^N. This was done by inserting a plasmid 
which coded for the production of the protein fragment into an £. coli strain and 
growing the genelically-modined bacterial strain in a limiting culture medium 
enriched with ^^NH4C1 and ^^C-glucose. 
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The isotopically enriched protein fragment was isolated from the culture 
medium, purified, and subsequently used as the basis for evaluating the binding of 
lest compounds. The procedures for these processes are described below. 

Human skin fibroblasts (ATCC No, CRL 1507) were grown and induced 

5 using the procedure described by Clark et al., Archiv. Biochem. and Biophvs . , 
24 1: 36-45 (1985). Total RNA was isolated from 1 g of cells using a Promega 
RNAgents® Total RNA Isolation System Kit (Cat.# Z51 10, Promega Corp., 2800 
Woods Hollow Road, Madison, WI 5371 1-5399) following the manufacturer's 
instructions. A 1 }4g portion of the RNA was heat-denatured at BO'C for five 

10 minutes and then subjected to reverse transcriptase PCR using a Gene Amp® RNA 
PCR kit (Cat.# N808-0017, Applied Biosystems/Perkin-Elmer. 761 Main Avenue, 
Norwalk, CT 06859-0156) following the manufacturer's instructions. 

Nested PCR was performed using first primers (A) GAAATGAAGAGTC 
TTCAA (SEQ ID NO:3) and (B) GCGTCCCAGGTTCTGGAG (SEQ ID NO:4) 

1 5 and thirty-five cycles of 94*C, two minutes; 45"C, two minutes; and 72*C three 
minutes. This was followed by reamplification with internal primers (C) 
ATACCATGGCCTATCCAT TGGATGGAGC (SEQ ID Na5) and (D) 
ATAGGATCCTTAGGTCTCAGGGGA GTCAGG (SEQ ID NO:6) using thirty 
cycles under the same conditions described immediately above to generate a DNA 

20 coding for amino acid residues 1-256 of human stromelysin. 

The PCR fragment was then cloned into PCR cloning vector pT7Blue(R) 
(Novagen, Inc.. 597 Science Drive, Madison, WI 5371 1) according to the 
manufacturer's instructions. The resulting plasmid was cut with Ncol and BamHI 
and the stromelysin fragment was subcloned into the Novagen expression vector 

25 pET3d (Novagen. Inc.. 597 Science Drive, Madison. WI 5371 1), again using the 
manufacturer's instructions. 

A mature stromelysin expression construct coding for amino acid residues 
81-256 plus an initiating methionine was generated from the 1-256 expression 
construct by PCR amplification. The resulting PCR fragment was first cloned into 

30 the Novagen pT7Blue(R) vector and then subcloned into the Novagen pET3d 
vector, using the manufacturer's instructions in the manner described above, to 
produce plasmid (pETST-83-256). This final plasmid is identical to that described 
by Qi-Zhuang et al.. Biochemistry . 3 1: 11231-11235(1992) with the exception 
that the present codes for a peptide sequence beginning two amino acids earlier, at 

35 position 81 in the sequence of human stromelysin. 



wo 97/18471 



18 



PCT/US96/18270 



Plasmid pErST-83-256 was transformed into £. coli strain 
BL21(DE3)/pLysS (Novagen, Inc., 597 Science Drive. Madison, WI 53711) in 
accordance with the manufacturer's instructions to generate an expression strain, 
BUZ l(DE3)/pLysS/pETST-255- 1 . 
5 A preculture medium was prepared by dissolving 1.698 g of 

Na2HP4^7H20. 0.45 g of KH2PO4. 0.075 g NaCl, 0.150 g ^^NH4C1. 0.300 
13C-glucose, 300 jiL of IM aqueous MgS04 solution and 15 /^L of aqueous 
CaCl2 solution in 150 mL of deionized water. 

The resulting solution of preculture medium was sterilized and transferred 
0 to a sterile 500 mL baffle flask. Immediately prior to inoculation of the preculture 
medium with the bacterial strain, 150 /^L of a solution containing 34 mg/mL of 
chloramphenicol in 100% ethanol and 1.5 mL of a solution containing 20 mg/mL 
of ampicillin were added to the flask contents. 

The flask contents were then inoculated with 1 mL of glycerol stock of 
15 genetically-modified £, Coli , strain BL21(DE3)/pLysS/pETST.255- 1. The flask 
contents were shaken (225 rpm) al 37'C until an optical density of 0.65 was 
observed, 

A fermentation nuuient medium was prepared by dissolving 1 13.28 g of 
Na2HP4«7H20, 30 g of KH2PO4. 5 g NaCI and 10 mL of 1% DF-60 antifoam 
20 agent in 9604 mL of deionized water. This solution was placed in a New 

Brunswick Scientific Miaos Fermenter (Edison, NJ) and sterilized at 12rc for 40 
minutes. 

Immediately prior to inoculation of the fermentation medium, the following 
pre-sterilized components were added to the fermentation vessel contents: 100 mL 
25 of a 10% aqueous solution of ^^NH4C1, 100 mL of a 10% aqueous solution of 
^^C-glucose, 20 mL of an aqueous IM solution of MgS04, 1 mL of an aqueous 
IM CaCl2 solution, 5 mL of an aqueous solution of thiamin hydrochloride ( 10 
mg/mL), 10 mL of a solution containing 34 mg/mL of chloramphenicol in 100% 
ethano! and L9 g of ampicillin dissolved in the chloramphenicol solution. The pH 
30 of the resulting solution was adjusted to pH 7.00 by the addition of an aqueous 
solution of4N H2SO4. 

The preculture of £. Coli , strain BL21(DE3)/pLysS/pETST-255-L from 
the shake-flask scale pnx^edure described above was added to the fermentor 
contents and cell growth was allowed to proceed until an optical density of 0.48 
35 was achieved During this process, the fermenter contents were automatically 
maintained at pH 7.0 by the addition of 4N H2SO4 or 
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4N KOH as needed. The dissolved oxygen content of the fermenter contents was 
maintained above 55% air saturation through a cascaded loop which increased 
agitation speed when the dissolved oxygen content dropped below 55%. Air was 
fed to the fermenter contents at 7 standard liters per minute (SLPM) and the culture 

5 temperature was maintained at 37'C throughout the process. 

The cells were harvested by centrifugation at 17,000 x g for 10 minutes at 
4'C and the resulting cell pellets were collected and stored at -85"C. The wet cell 
yield was 3.5 g/L. Analysis of the soluble and insoluble fractions of cell lysates by 
sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) revealed 

10 that approximately 50% of the ^^N-stromelysin was found in the soluble phase. 

The isotopically-labeled stromelysin fragment prepared as described above 
was purified employing a modification of the technique described by Yeetal . , 
Biochemistry, 31: 11231-11235(1992). 

The harvested cells were suspended in 20 mM Tris-HCl buffer (pH 8.0) 

1 5 sodium azide solution containing 1 mM MgCl2> 0.5 mM ZnCl2. 25 units/mL of 
Benzonase® enzyme, and an inhibitor mixture made up of 4-(2-aminoethyl)- 
benzenesulfonyl fluoride ("AEBSF'), Leupeptin®. Aprotinin®, and Pepstatin® 
(all at concentrations of 1 /<g/mL. AEBSF, Leupeptin®, Aprotinin®, and 
Pepstatin® are available from American International Chemical. 17 Strathmore 

20 Road, Natick, MA 01760.) 

The resulting mixture was gently stirred for one hour and then cooled to 
4*0 The cells were then sonically disrupted using a 50% duty cycle. The 
resulting lysate was centrifuged at 14,000 rpm for 30 minutes and the pellet of 
insoluble fraction frozen at -80'C for subsequent processing (see below). 

25 Solid sunmonium sulfate was added to the supernatant to the point of 20% 

of saturation and the resulting solution loaded onto a 700 mL phenyl sepharose fast 
flow ("Q-Sepharose FF*) column (Pharmacia Biotech.. 800 Centennial Ave.. P. 
O. Box 1327, Piscataway. NJ 08855). Prior to loading, the sepharose column 
was equilibrated with 50 mM Tris-HCl buffer (pH 7.6 at 4''G). 5 mM CaCl2, and 

30 1 M (NH4)2S04. The loaded column was eluted with a linear gradient of 
decreasing concentrations of aqueous (NH4)2S04 (from 1 down to 0 M) and 
increasing concentrations of aqueous CaCl2 (f-rom 5 to 20 mM) in Tris-HCl buffer 
at pH 7.6. 



35 
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The active fractions of eluate were collected and concentrated in an Amicon 
stirred cell (Amicon, Inc., 72 Cherry Hill Drive, Beverly, MA 01915). The 
concentrated sample was dialy7:ed overnight in the starting buffer used y/'\ih the Q- 
Scpharose FF column, 50 mM Tris-HCl (pH 8.2 at 4'C) with 10 mM CaCl2. 
5 The dialyzed sample was then loaded on the Q-Sepharose FF column and 

eluted with a linear gradient comprising the starting buffer and 200 mM NaCl. The 
purified soluble fraction of the isotopically-labeled stromelysin fragment was 
concentrated and stored at 4°C. 

The pellet was solubilized in 8M guanidine-HCL The solution was 
1 0 centrifuged for 20 minutes at 20,000 rpm and the supernatant was added dropwise 
to a folding buffer comprising 50 mM Tris-HCl (pH 7.6). 10 mM CaCl2 0.5 mM 
ZnCl2 and the inhibitor cocktail of AEBSF. Leupeptin®. Aprotinin®. and 
Pepstatin® (all at concentrations of 1 i^glmL). The volume of folding buffer was 
ten times that of the supernatant. The mixture of supernatant and folding buffer 
1 5 was centrifuged at 20,000 rpm for 30 minutes. 

The supernatant from this centrifugation was stored at4*C and the pellet 
was subjected twice to the steps described above of solubilization in guanidine- 
HCl, refolding in buffer, and centrifugation. The final supematants from each of 
the three centrifugations were combined and solid ammonium sulfate was added to 
20 the point of 20% saturation. The resulting solution thus derived from the insoluble 
fraction was subjected to purification on phenyl Sepharosc and Q-Sepharose as 
described above for the soluble fraction. 

The purified soluble and insoluble fractions were combined to produce 
about 1.8 mg of purified isotopically-labeled stromelysin 81-256 fragment per 
25 gram of original cell paste. 

B. Human papillomavirus (HPV) E2 Inhibitors 

The papillomaviruses are a family of small DNA viruses that cause genital 

warts and cervical carcinomas. The E2 protein of HPV regulates viral transcription 

30 and is required for viral replication. Thus, molecules that block the binding of E2 

to DNA may be useful therapeutic agents against HPV. The protein rather than the 

DNA was chosen as a target, because it is expected that agents with greater 

selectivity would be found that bind to the protein rather than the DNA. 

The DNA-binding domain of human papillomavirus E2 was cloned from 

35 the full length DNA that codes for E2 using PCR and overexpressed in bacteria 

15 

using the T7 expression system. Uniformly N-labeied protein was isolated from 
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bacteria grown on a minimal medium containing N-labclcd protein was isolated 

from bacteria grown on a minimal medium containing ^^N-labeied ammonium 

chloride. The protein was purified from the bacterial cell lysate using an S- 

sepharose FastRow column pre-equilibrated with buffer (50 mM Tris, 100 mM 

NaCl, 1 mM EDTA, pH = 83), 

The protein was eluted with a linear gradient of 100-500 mM NaCl in 

buffer, pooled, and applied to a Mono-S column at a pH = 7.0. The protein was 

eluted with a salt gradient (100-500 mM), concentrated to 0.3 mM, and exchanged 

into a TRIS (50 mM, pH = 7.0 buffered H2CyD20 (9/1) solution containing 

sodium azide (0.5%). 

C. EAE 

Uniformly ^^N-labeled Ras-binding domain of the RAF protein was 
prepared as described in Emerson et al.. Biochemistry . 34 (21): 691 1-6918 
(1995). 

D. FKBP 

Uniformly ^^N-labeled recombinant human FK binding protein (FKBP) 
was prepared as described in Logan et al., J. Mol. Biol. . 236: 637-648 (1994). 

Example 2 

Screening Compounds Using Two-Dimensional — N/- H NMR Correlation 
Spectral Analysis 

The catalytic domain of stromelysin was prepared in accordance with the 
procedures of Example 1. The protein solutions used in the screening assay 
contained the uniformly ^^N-labeled catalytic domain of stromelysin (0.3 mM), 
acetohydroxamic acid (500 mM). CaCl2 (20 mM), and sodium azide (0.5%) in a 
H2O/D2O (9/1) TRIS buffered solution (50 mM. pH=7.0), 

Two-dimensional ^^N/^H NMR spectra were generated at 29*C on a 
Bmker AMX500 NMR spectrometer equipped with a triple resonance probe and 
Bruker sample changer. The ^^N/^H HSQC spectra were acquired as 80 x 1024 
complex points using sweep widths of 2000 Hz (^^N. t^) and 8333 Hz (^H, t2). 
A delay of 1 second between scans and 8 scans per free induction dccay(fid) were 
employed in the data collection. All NMR spectra were processed and analyzed on 
Silicon Graphics computers using in-house- written software. 
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A first two-dimensional N/ H NMR correlation speclnim was acquired 
for the ^^N-labeled stromelysin target molecule as described above. The 
stromelysin target was then exposed to a database of test compounds. Stock 
solutions of the compounds were made at 100' mM and 1 M. In addition, a 
5 combination library was prepared that contained 8-10 compounds per sample at a 
concentration of 100 mM for each compound. 

The pH of the 1 M stock solution was adjusted with acetic acid and 
ethanolamine so that no pH change was observed upon a 1/10 dilution with a 100 
mM phosphate buffered solution (pH = 7.0). It is important to adjust the pH, 
10 because small changes in pH can alter the chemical shifts of the biomolecules and 
complicate the interpretation of the NMR data. 

The compounds in the database were selected on the basis of size 
(molecular weight ~ 100-300) and molecular diversity. The molecules in the 
collection had different shapes (e.g.. flat aromatic rings(s), puckered aliphatic 
1 5 rings(s), straight and branched chain aliphatics with single, double, or triple 
bonds) and diverse functional groups (e.g.. carboxylic acids, esters, ethers, 
amines, aldehydes, ketones, and various heterocyclic rings) for maximizing the 
possibility of discovering compound that interact with widely diverse binding sites. 
The NMR samples were prepared by adding 4 lA of the DMSO slock 
20 solution of the compound mixtures that contained each compound at a 

concentration of 100 mM to 0.4 ml H2CyD20 (9/1) buffered solution of the 
uniformly ^^N-labeled protein. The final concentration of each of the compounds 
in the NMR sample was about 1 mM. 

In an initial screen, two compounds were found that bind to the catalytic 
25 domain of stromelysin. Both of these compounds contain a biaryl moiety. Based 
on these initial hits, structurally similar compounds were tested against 
stromelysin. The structure of those biaryl compounds is represented by the 
structure I . below. (See Table 1 for definitions of R1-R3 and A 1-A3) . 

30 "2 

In the second round of screening, binding was assayed both in the absence 
and in the presence of saturating amounts of acetohydroxamic acid (500 mM). 
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Many of the biaryl compounds were found to bind the catalytic domain of 
slromelysin. FIG. 4 shows a representative two-dimensional ^^N/^H NMR 
correlation spectrum before and after exposure of stromelysin to a biaryl test 
compound. It can be seen from FIG, 4 that the compound caused chemical shifts 

5 of ^^N-sites such as those designated W124. T187, A 199 and 0204, 

These sites correspond to a tryptophan (Trp) residue at position 124. a 
threonine (Thr) at position 187, an alanine (Ala) at positionl99, and a glycine 
(Gly) at position 204 of SEQ ID NO. 1. FIG. 9 shows the correlation between the 
NMR binding data and a view of the NMR-derived three-dimensional structure of 

10 the catalytic domain of stromelysin. The ability to locate the specific binding site of 
a particular ligand is an advantage of the present invention. 

Some compounds only bound to stromelysin in the presence of hydroxamic 
acid. Thus, the binding affinity of some compounds was enhanced in the presence 
of the hydroxamic acid (i. e. cooperative). These results exemplify another 

15 important capability of the present screening assay: the ability to identify 
compounds that bind to the protein in the presence of other molecules. 

Various biaryl compounds of structure I were tested for binding to 
stromelysin at differing concentrations. The ^^N/^H spectra generated at each 
concentration were evaluated to quantify differences in the spectra as a function of 

20 compound concentration. A binding or dissociation constant (KD)was calculated, 
using standard procedures well known in the art, from those differences. The 
results of this study are shown in Table 1. The values for R1-R3 and A 1-A3 in 
Table 1 refer to the corresponding positions in the structure 1, above. 

25 Table 1 



Compound 
No. 


Rl 


R2 


R3 


Ai 


A2 


A3 


Kd(itiM) 
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OH 
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C 
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CH2OH 
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C 
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3.2 
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Br 
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OH 
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1.3 
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" N 


N 
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1.6 
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CHO 


H 
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c 


c 


1.7 


6 


OCH3 


NH2 


H 


C 


c 


c 


0.4 


7 


H 


H 


H 


N 


c 


c 


0.2 
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8 


OCOCH3 


H 


H 


C 


c 


c 


0.3 


9 


OH 


H 


OH 


C 


c 


c 


0.16 


10 


H 


H 


H 


N 


c 


N 


0.4 


11 


OH 


H 


H 


C 


c 


c 


0.3 


12 


OH 


H 


CN 


C 


c 


c 


0.02 



The data in Table 1 show the utility of a process of the present invention in 
determining dissociation or binding constants between a ligand and a target 
molecule. 

5 Another advantage of an NN4R screening assay of the present invention is 

15 1 

the ability to correlate observed chemical shifts from the two-dimensional N/ H 
NMR correlation spectra with other spectra or projections of target molecule 
configuration. The results of a representative such correlation are shown in FIG. 
9. which depicts regions within the polypeptide at which binding with the substrate 
10 molecule is most likely occurring. In this Figure, the apparent binding regions in 
stromelysin are shown for Compound I (from Table 1), 

Compounds from the database were screened in a similar manner for 
binding to the DNA-binding domain of the E2 protein. Those compounds had the 
structure II below, where R1-R4 and A arc defined in Table 2. 




NMR experiments were performed at 29'C on a Brukef AMXSOO NMR 
spectrometer equipped with a triple resonance probe and Bruker sample changer. 
20 The ^^N-/^H HSQC spectra were acquired as 80 x 1024 complex points using 
sweep widths of2000 Hz (^^N,ti) and 8333 Hz(^H, t2). Adelay of 1 second 
between scans and 4 scans per free induction decay were employed in the data 
collection. All NMR spectra were processed and analyzed on Silicon Graphics 
computers. 

25 FIGs. 2 and 3 show representative two-dimensional ^^N/^H NMR 

correlation spectra before and after exposure of the DNA-binding domain of E2 to 
a first and second test compound, respectively. 
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It can be seen from FIG. 2 that the first test compound caused chemical 
shifts of ^^N-sites such as those designated 115, Y21, R22 and L23. Those sites 
correspond to an isoleucine (He) residue at position 15, a tyrosine residue (Tyr) at 
position 21, an arginine (Arg residue at position 22 and a leucine (Leu) residue at 
5 position 23 of SEQ ID NO. 6. 

It can be seen from RG. 3 that the second test compound caused chemical 
shifts in the particular ^^N-sites designated 16, Gil, H38. and T52. Those sites 
correspond to an isoleucine (He) residue at position 6, a glycine (Gly) residue at 
position 11, a histidine (His) residue at position 38 and a threonine (Thr) at 
10 position 52 of SEQ ID NO. 6. 

HGs. 7 and 8 show the correlation between those NMR binding data and a 
view of the NMR-deri ved three-dimensional structure of the DNA-binding domain 
of E2. 

Several structurally similar compounds caused chemical shift changes of 
15 the protein signals when screened at a concentration of 1 mM. Two distinct sets of 
amide resonances were found to change upon the addition of the compounds: one 
set of signals corresponding to amides located in the B-barrel formed between the 
two monomers and a second set corresponding to amides located near the DNA- 
binding site. 

20 For example, compounds containing two phenyl rings with a carboxylic 

acid attached to the carbon linking the two rings only caused chemical shift changes 
to the amides in the DNA-binding site. In contrast, benzophenones and 
phenoxyphenyl-containing compounds only bound to the B-barrel. Other 
compounds caused chemical shift changes of both sets of signals but shifted the 

25 signals in each set by different amounts, suggesting the presence of two distinct 
binding sites. 

By monitoring the chemical shift changes as a function of ligand 
concentration, binding constants for the two binding sites were also measured. 
The results of those studies are summarized below in Table 2. 
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Table 2 



Comp. 
No. 


A 




R2 


R3 


R4 


DNA 


P-band 
KD(niM) 


Filler 
binding 
assay 


13 


CO 
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H 


OH 


>50 


0.6 




14 
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H 


CH2OH 


>50 


2.0 
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,a 
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COO 
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2.0 


>50 


-f 


16 


_a 


CI 


CI 


COO 
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>50 


+ 


17 


_a 
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CH2COO 
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4.2 


4.9 


-I- 


18 


.a 
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CH=CHCOO 


H 


1-2 


6.2 


+ 


19 
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CH2CH2CH(CH3) 
-CH2COO 
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0.5 


0.2 




20 


0 


H 


H 


COCH2CH2COO 


H 


2.7 


4.8 


+ 



^ a dash (-) for A indicates no atom (i.e. byphenyl linkage) 



5 Uniformly N-labeled Ras-binding domain of the RAF protein was 

prepared as described in Example 1 and screened using two-dimensional N/ H 
NMR correlation spectral analysis in accordance with the NMR procedures 
described above. The results of a representative study are shown in FIG. 5, which 
depicts two-dimensional ^^N/^H NMR correlation spectra both before and after 

10 exposure to a test compound. 

Uniformly ^^N-labeled FKBP was prepared as described in Example 1 and 
screened using two-dimensional ^^N/^H NMR correlation spectral analysis in 
accordance with the NMR procedures described above. The results of a ^ 
representative study are shown in FIG. 6, which depicts two-dimensional N/ H 

15 NMR correlation spectra both before and after exposure to a test compound. 

Example 3 

rom parison of NMR. FnTvmatic. Filte r Rinding and Gel Shift Screening Assays 
Studies were performed to compare binding constants of ligands to various 
20 biomolecules. detennincd by the NMR method of the present invention, to similar 
results obtained from prior art methods. 

In a first study, binding constants were determined, both by the NMR 
method of the present invention, and by a prior art enzymatic assay. The target 
molecule was the catalytic domain of siromelysin prepared in accordance with the 
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procedures of Example L The NMR binding constants. Kd, were derived using 
two-dimensional ^^N/^H NMR correlation spectroscopy as described in Example 
2. The Kd values so obtained were compared to an inhibition constant Ki as 
determined in an enzymatic assay. 

5 The enzymatic assay measured the rate of cleavage of a fluorogenic 

substrate by following the fluorescence increase upon peptide cleavage which 
causes a separation between the fluorophore and quencher. Enzymatic activity was 
measured using a matrix of different concentrations of acetohydroxamic acid and 
biaryl compounds. The assay is a modification of the method described by H. 

10 Weingarten, et al in Anal. Biochcm. , 147: 437-440 (1985) employing the 

fluorogenic substrate properties described by E. Matayoshi, et al in Science : 247: 
954-958(1990). 

Eight acetohydroxamic acid concentrations were used ranging from 0.0 to 
1 .0 M, and six compound concentrations were used, resulting in a total of 48 
15 points. Individual compound concentration varied due to solubility and potency. 

All NMR measurements were performed in the presence of 500 mM 
acetohydroxamic acid, except for the titration of acetohydroxamic acid itself. 
Dissociation constants were obtained from the dependence of the observed 
chemical shift changes upon added ligand. Inhibition constants were then obtained 
20 from the inhibition data using standard procedures. 

The results of these studies are summarized below in Table 3, which shows 
the comparison of NMR-derived dissociation constants (Kd) with inhibition 
constants measured in the enzyme assay (Kj) using a fluorogenic substrate. 

25 Table 3 



Compound 
No. 


NMR Kd 
(mM) 


Assay Kl 
(mM) 


4 


1.6 


7.4 


7 


0.17 


0.32 


9 


0.16 


0.70 


10 


•0.40 


1.8 


12 


0.02 


0.11 


Acetohydroxamic acid 


17.0 


21.1 
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The data in Table 3 show that a NMR process of the present invention 
provides a rapid, efficient and accurate way of determining dissociation or binding 
constants of ligands to target biomolecules. Comparison of the binding constants 
determined by the two methods result in the same ranking of potencies of the 
5 compounds tested. That is. while the values for a given substrate as determined by 
the two methods are not equal, Ihey are proportional to one another. 

In a second study, the results for binding of the DNA-binding domain of 
E2 to its target DNA were obtained by prior art methods and compared with results 
obtained by the method of the present invention. The target was the DNA-binding 
10 domain of E2, prepared in accordance with the procedures of Example 1 . NMR 
screening assays and NMR processes for determining ligand dissociation constants 
were performed as set forth above in Example 2. 

The binding constant from the NMR process was compared to the results 
of a physical, filter binding assay that measured binding of DNA to the target. The 
15 high-throughput filter binding assay was performed using E2. prepared according 
to Example 2 above. The ^'^P-labeled DNA construct comprised a 10329 base 
pair plasmid formed by inserting the HPV-1 1 genome, containing three high 
affinity and one low affinity E2 binding sites, into the PSP-65 plasmid (Promega. 
Madison, Wl), 

20 The binding affinities at the different sites as detennined by NMR were 

compared for a subset of the compounds to the inhibition of E2 binding to DNA as 
measured in the filter binding assay. As shown in Table 2 above, the activities 
determined in the filter binding assay correlated closely with the binding affinities 
calculated from the amides of the DNA-binding site but not to the affinities 

25 measured for the B-barrel site. This is consistent with the relative locations of each 
site. 

In an alternative study, a comparison of the NMR-determined binding 
results was made with similar results obtained by a prior art gel-shift assay using 
techniques well known in the art. The gel-shift assay was performed using a GST 

30 fusion protein which contained full length E2 and a 33p-iabeled 62 base pair DNA 
fragment containing two E2 binding sites. 

The method identified numerous compounds which gave positive results in 
the gel-shift assay. Some of these positive results, however, were believed to be 
due to binding to the DNA, since in these cases, no binding to the E2 protein was 

35 observed using the NMR method of this invention. These compounds were 
shown to indeed bind to DNA rather than to E2, as evidenced by changes in the 
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chemical shifts of the DNA rather than the protein upon the addition of the 
compounds. These data show that yet another advantage of the present invention 
is the ability to minimize the occurrence of false positives. 

5 Example 4 

Design of a potent, non-peptide inhibitor of stromelvsin 

Studies were performed to design new ligands that bound to the catalytic 
domain of stromelysin. Because stromelysin undergoes autolysis, an inhibitor was 
sought to block the degradation of stromelysin. That inhibitor would facilitate the 
10 screening of other potential ligands that bind to other sites on the enzyme. 

The criteria used in selecting compounds in the screening for other binding 
sites was based primarily on the size of the ligand. The smallest iigand was sought 
that had enough solubility to saturate (>98% occupancy of enzyme) and inhibit the 
enzyme. 

1 5 The cloning, expression, and purification of the catalytic domain of 

stromelysin was accomplished using the procedures set forth in Example 1. An 
initial step in the design of the new ligand was the identification of a first ligand 
that bound to the stromelysin target Such identification was carried out in 
accordance with a two-dimensional ^^N/^H NMR correlation screening process as 

20 disclosed above. 

A variety of hydroxamic acids of the general formula R-(CO)NHOH were 
screened for binding to stromelysin using the procedures set forth in Example 2. 
Of the compounds tested, acetohydroxamic acid lCH3(CO)NHOH] best satisfied 
the selection criteria: it had a binding affinity for stromelysin of 17 mM and had 

25 good water solubility. At a concentration of 500 mM, acetohydroxamic acid 

inhibited the degradation of the enzyme, allowing the screening of other potential 
ligands. 

The second step in the design process was the identification of a second 
ligand that bound to the target stromelysin at a site different from the binding site of 
30 acetohydroxamic acid. This was accomplished by screening compounds for their 
ability to bind stromelysin in the presence of saturating amounts of 
acetohydroxamic acid. Details of procedures' and results of this second 
identification step are set forth above in Example 2, 

The compound identified as a second ligand from these studies and used in 
35 subsequent design steps was the compound designated as Compound #4 in Table 1 
(See Example 2). 
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The next step in the design process was to constnict a ternary complex of 
the target stromelysin, the first hgand and the second ligand. This was 
accomplished by exposing the stromelysin target to the two ligands under 
conditions that resulted in complex formation. The three-dimensional structure of 
5 the ternary complex was then determined using NN4R spectroscopy as described 
below. 

The ^H, ^^C, and ^^N backbone resonances of stromelysin in the ternary 
complex were assigned from an analysis of several 3D double- and triple- 
resonance NMR spectra (A. Bax. ef fl/.. Acc. Chem. Res. . 26: 131-138(1993)). 
10 The resonances of adjacent spin systems were identified from an analysis of 
three-dimensional (3D) HNCA (L. Kay etal., J. Magn. Reson., 89: 496-514 
(1990)) and HN(CO)CA (A. Bax, etaL, J. Bio. NMR , 1 : 99 (1991)) spectra 
recorded with identical spectral widths of 1773 Hz (35.0 ppm). 3788 Hz (30. 1 
ppm). and 8333 Hz (16.67 ppm) in the Fi(^^N), F2('^C) and F3(^H) 

15 di mensions . respec ti vely . 

The data matrix was 38(t i) x 48(t2) x 1024{t3) complex points for the 
HNCA spectrum, and 32(1 1) x 40(t2) x 1024{t3) complex points for the 
HN(CO)CA spectrum. Both spectra were acquired with 16 scans per increment. A 
3D CBCA(CO)NH spectrum (S. Grzesiek, etai, J. Am. Chem. Soc., 114: 6261- 
20 6293 (1992)) was collected with 32(ti, ^^N) x 48(t2. ^^C) x 1024(t3. ^H) 

complex points and 32 scans per increment. Spectral widths were 1773 Hz (35.0 
ppm), 7575.8 Hz (60.2 ppm). and 8333 Hz (16.67 ppm) in the ^^N, ^ C and H 
dimensions, respectively. 

For all three spectra, the ^H carrier frequency was set on the water 
^ resonance and the ^^N carrier frequency was at 1 19. 1 ppm. The ^^C carrier 

frequency was set to 55.0 ppm in HNCA and HN(CO)CA experiments, and 46.0 
ppm in the CBCA(CO)NH experiment. 

The backbone assignments were confirmed from an analysis of the 
crosspeaks observed in an ^^N-separated 3D NOESY-HSQC spectrum and a 3D 
30 HNHA-J spectrum. The ^^N-separated 3D NOESY-HSQC spectrum (S. Fesik. 
etai^ J. Magn. Reson. , 87: 588-593 (1988)); D. Marion, etai, J. Am. Chem. 
Soc. Ill: 1515-1517 (1989)) was collected with a mixing time of 80 ms. A total 
of 68(ti. ^^N) X 96(t2, ^H) x 1024(13, ^H) complex points with 16 scans per 
increment were collected, and the spectral widths were 1773 Hz (35.0 ppm) for the 
35 ^^N dimension, 6666.6 Hz (12, ^H, 13.3 ppm), and 8333 Hz (16.7 ppm) for the 
^H dimension. 
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The 3D HNHA-J spectrum (G. Vuister, et al, J. Am. Chem. Soc. . 1 1 5: 

7772-7777 (1993)), which was also used to obtain ^JHNHa coupling constants, 

was acquired with 35(ti, ^^N)x64{l2, *H)xl024(t3, ^H) complex points and 

32 scans per increment. Spectral widths and carrier frequencies were identical to 

5 those of the ^^N-separated NOESY-HSQC Sfjectrum, Several of the h'^ signals 

were assigned using the HNHB experiment. The sweep widths were the same as 

in the ^^N-separated NOESY-HSQC spectrum that was acquired with 32(ti. ^^N) 

X 96(t2, ^H) X 1024 (13, ^H) complex points. 
1 13 

The H and C chemical shifts were assigned for nearly all sidechain 

10 resonances. A 3D HCCH-TOCSY spectrum (L. Kay, et a/.. J, Magn. Reson. . 

101b: 333-337 ( 1993)) was acquired with a mixing time of 13 ms using the 

DIPSI-2 sequence (S. Rucker, et al, Mol. Phvs. . 68: 509 (1989)) for ^"^C 

isotropic mixing. A total of 96 (ti, ^^C) x96(t2, ^H)x 1024{t3. ^H) complex 

data points were collected with 16 scans per increment using a spectral width of 

15 10638 Hz (70.8 ppm. wi). 4000 Hz (6.67 ppm. W2). and 4844 (8.07 ppm. W3). 

13 

Carrier positions were 40 ppm, 2.5 ppm, and at the water frequency for the C , 
indirectly detected ^H, and observed dimensions, respectively. 

13 

Another 3D HCCH-TOCSY study was performed with the C carrier at 
122.5 ppm to assign the aromatic residues. The spectra were collected with 
20 36(ti, C) x48(t2,'^H) x 1024 (13. ^H) complex points with specU-al widths of 
5263 Hz (35.0 ppm. wi), 3180 Hz (5.30 ppm, W2). and 10.000 (16.7 ppm. W3). 

Carrier positions were 122.5 ppm, 7.5 ppm. and at the water frequency for the 

13 11 
C, indirectly detected H. and observed H dimensions, respectively. 

A ^^C-separated 3D NOESY-HMQC spectrum (S. Fesik, et ai, J. Magn. 

25 Reson. . 8 7: 588-593 (1988)); D. Marion, etai, J. Am. Chem. Soc. . 1 1 1: 1515- 

13 

1517(1989)) was recorded using a mixing time of 75 ms. A total of 80(ti, C) 

X 72 (t2, ^H) X 1024 (13, ^H) complex data points with 16 scans per increment 

were collected over spectral widths of 10638 Hz (70.49 ppm. w 1), 6666.6 Hz 

(13.3 ppm, W2). and 8333.3 Hz (16.67 ppm. W3). The ^H carrier frequencies 

13 

30 were set to the water resonance, and the C carrier frequency was placed at 40.0 
ppm. 

Stereospecific assignments of methyl groups of the valine and leucine 

residues were obtained by using a biosynthetic approach (Neri et ai, Biochem. . 

28: 7510-7516 (1989)) on the basis of the ^'^C-^'^C one-bond couphng pattern 

1 13 

35 observed in a high-resolution H, C-HSQC spectrum (G. Bodenhausen, et ai, 
J. Chem. Phvs. Lett. , 69: 185-189 (1980)) of a fractionally ^-^C-labeled protein 
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sample. The spectrum was acquired with 200(*''c. ti) x2048('H, 12) complex 

points over spectral widths of 5000 Hz (39.8 ppm, ^^C) and 8333 Hz (16.7 ppm, 

^H). Carrier positions were 20.0 ppm for the ^^C dimension, and at the water 

frequency for the dimension. or> 12^ 

5 To detect NOEs between the two ligands and the protein, a 3D 

filtered. ^''c-edited NOESY spectrum was collected. The pulse scheme consisted 

of a double ^^C-filter sequence (A. Gemmeker, et al, J. Magn. Reson.. 9 6: 199- 

204 (1992)) concatenated with a NOESY -HMQC sequence (S. Fesik, et al.. L 

Mapn. Reson. . 87: 588-593 (1988)); D. Marion, etai. J. Am. Chem. Soc.. 1 1 1: 

10 1515-1517 (1989)) . The spectrum was recorded with a mixing time of 80 ms, 
and a total of 80 (ti, ^^C) x 80 (t2. ^H) x 1024 (t3. ^H) complex points with 16 
scans per increment. Spectral widths were 8865 Hz ( 17.73 ppm. w i), 6667 Hz 
( 13.33 ppm, W2). and 8333 Hz (16.67 ppm, W3), and the carrier positions were 
40.0 ppm for the carbon dimension and at the water frequency for both proton 

15 dimensions. 

To identify amide groups that exchanged slowly with the solvent, a series 
of ^H. ^^N-HSQC spectra (G. Bodcnhausen. et ai, I Cj^hem. Phvs. Lett.. 69: 
185-189 (1980)) were recorded at 25'C at 2 hr intervals after the protein was 
exchanged into D2O. The acquisition of the first HSQC spectrum was started 2 
20 hrs. after the addition of D2O. 

All NMR spectra were recorded at 25'C on a Bruker AMX500 or AMX600 
NMR spectrometer. The NMR data were processed and analyzed on Silicon 
Graphics computers. In all NMR experiments, pulsed field gradients were applied 
where appropriate as described (A. Bax. et ai, J, Magn. Reson.. 9 9: 638 ( 1992)) 
25 to afford the suppression of the solvent signal and spectral artifacts. Quadrature 
detection in indirecUy detected dimensions was accomplished by using the Stales- 
TPPI method (D. Marion, et ai, J Am. Chem. Soc.. 1 1 1: 1515-1517 ( 1989)). 
Linear prediction was employed as described (E. Olejniczak, et ai, J. Magn. 
Reson.. 87: 628-632 (1990)). 
30 The derived three-dimensional structure of the ternary complex was then 

used to define the spatial orientation of the first and second li gands to each other as 
well as to the target stromelysin molecule. 

Distance restraints derived from the NOE data were classified into six 
categories based on the NOE cross peak intensity and given a lower bound of 1.8 
35 A and upper bounds of 2.5 A, 3.0 A, 3.5 A. 4.0 A, 4.5 A, and 5.0 A. 

respectively. Restraints for <|> torsional angles were derived from JHNHa 
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coupling constants measured from the 3D HNHA-J spectrum (G. Vuister, et qL, L. 
Am. Chem. Soc. > 1 1 5: llll'llll (1993)). The ^ angle was restrained to 
120%±40% for ^JHNHa > 8.5 Hz, and 60%±40% for ^JHNHa < 5 Hz. 

Hydrogen bonds, identified for slowly exchanging amides based on initial 

5 structures, were defined by two restraints: 1.8-2.5 A for the H-0 distance and 1.8- 
3,3 A for the N-O distance. Structures were calculated with the X-PLOR 3. 1 
program (A. BrUngcr. "XPLOR 3.1 Manual." Yale University Press, New Haven. 
1992) on Silicon Graphics computers using a hybrid distance geometry-simulated 
annealing approach (M. Nilges, et a/., FEBS Lett. . 229: 317-324 (1988)). 

,0 A total of 1032 approximate interproton distance restraints were derived 

from the NOE data. In addition, 21 unambiguous intermolecular distance restraints 
were derived from a 3D 12C-filtered, 13C-edited NOESY spectrum. Of the 1032 
NOE restraints involving the protein, 341 were intra-residue. 410 were sequential 
or short-range between residues separated in the primary sequence by less than five 

15 amino acids, and 281 were long-range involving residues separated by at least five 
residues. 

In addition to the NOE distance restraints, 14 <j> dihedral angle restraints 
were included in the suiicture calculations that were derived from three-bond 
coupling constants ( JHNHa) determined from an HNHA-J spectrum (G. 
20 Viioster, etd,, J. Am. Chem. Soc. , 1 1 5: 1H1-11TI ( 1993)). The experimental 
restraints also included 120 distance restraints corresponding to 60 hydrogen 
bonds. The amides involved in hydrogen bonds were identified based on their 
characteristically slow exchange rate, and the hydrogen bond partners from initial 
NMR structures calculated without the hydrogen bond restraints. The total number 
25 of non-redundant, experimentally-derived restraints was 1 166. 

The structures were in excellent agreement with the NMR experimental 
restraints. There were no distance violations greater than 0.4 A. and no dihedral 
angle violations greater than 5 degrees. In addition, the simulated energy for the 
van der Waals repulsion term was small, indicating that the structures were devoid 
30 of bad inter-atomic contacts. 

The NN4R structures also exhibited good covalent bond geometry, as 
indicated by small bond-length and bond-angle deviations from the corresponding 
idealized parameters. The average atomic root mean square deviation of the 8 
structures for residues 93-247 from the mean coordinates was 0.93 A for backbone 
35 atoms (C^ N, and C). and 1.43 A for all non-hydrogen atoms. 
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A ribbon plot of the ternary complex involving siromclysin, 
acelohydroxamic acid (the first ligand). and the second ligand is shou^n in Fig 10. 
The structure is very similar to the global fold of other matrix metal loproteinases 
and consists of a five-stranded B-sheet and three a-heliccs. 

5 The catalytic zinc was located in the binding pocket. It was coordinated to 

three hislidines and the two oxygen atom of acetohydroxamic acid, A biaryl group 
of the second ligand was located in the ST pocket between the second helix and 
the loop formed from residues 218-223. This deep and narrow pocket is lined 
with hydrophobic residues which make favorable contacts with the ligand. 

10 Based on the three-dimensional structure of the ternary complex as 

determined above and the structure/activity relationships observed for the binding 
to stromelysin of structural analogs of the second ligand (i.e., other biaryl 
compounds), new molecules were designed that linked together the 
acetohydroxamic acid to biaryls. 

15 As shown in Table 4 below, the initial biaryls chosen contained an oxygen 

linker and the absence or presence of CN para to the biaryl linkage. Initial linkers 
contained varying lengths of methylene units. Means for linking compounds with 
linkers having varying lengths of methylene units arc well known in the arl 



Table 4 
H 




Compound 


X 


R 


Stromelysin 
I Inhibition 


21 


(CH2)2 


H 


0.31 fiM 


22 


(CH2)3 


H 


110 }4M 


23 


(CH2)4 


H 


38%@100/<M 
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24 


(CH2)5 


H 


43%@100;<M 


25 


(CH2)2 


CN 


0.025 iaM 


26 


(CH2)3 


CN 


3.4 |<M 


27 


(CH2)4 


CN 


3.5 iaM 


28 


(CH2)5 


CN 


1.7 ixM 



As expected based on the better binding of the CN substituted biaryls to 
stromelysin, the CN derivatives exhibited better stromelysin inhibition. The 
compound that exhibited the best inhibition of stromelysin contained a Hnker with 
5 two methylene units. 

The present invention has been described with reference to preferred 
embodiments. Those embodiments are not limiting of the claims and specification 
in any way. One of ordinary skill in the art can readily envision changes, 
modifications and alterations to those embodiments that do not depart from the 
10 scope and spirit of the present invention. 
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SEQUENCE LISTING 



(1) GENERAL INFORMATION: 

5 

(i) APPLICANT: Fesik, Stephen W. 

Hajduk, Philip J. 

(ii) TITLE OF INVENTION: Use of Nuclear Magnetic 
10 Resonance to 

Identifyify Ligands to Target Biomolecules 
(iii) NUMBER OF SEQUENCES: 6 

15 (iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: Steven F- Weinstock, Dept. 377 

AP6D, Abbott 
Laboratories 

(B) STREET: 100 Abbott Park Road 
20 (C) CITY: Abbott Park 

(D) STATE: Illinois 

(E) COUNTRY: USA 
(P) ZIP: 60064-3500 

25 (V) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS /MS-DOS 

(D) SOFTWARE: PatentIn Release #1.0, Version 

30 #1.30 

(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER: 

(B) FILING DATE: 

35 (C) CLASSIFICATION: 

(viii) ATTORNEY /AGENT INFORMATION: 

(A) NAME; Janssen, Jerry F. 

(B) REGISTRATION NUMBER: 29,175 



40 



45 



<ix) TELECOMMUNICATION INFORMATION: 

(A) TELEPHONE: (708) 937-4558 

(B) TELEFAX: (708) 938-7742 



(2) INFORMATION FOR SEQ ID N0:1: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 174 amino acids 
50 (B) TYPE: amino acid 

(C) ' STRANDEDNESS : single 

(D) TOPOLOGY: linear 
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(ii) MOLECULE TYPE: peptide 

(xi) SEQUENCE DESCRIPTION: SEQ ID N0:1: 

5 

Phe Arg Thr Phe Pro Gly He Pro Lys Trp Arg Lys Thr 
1 5 10 

His Leu Thr Tyr Arg He Val Asn Tyr Thr Pro Asp Leu 
10 15 20 25 

Pro Lys Asp Ala Val Asp Ser Ala Val Glu Lys Ala Leu 
30 35 

15 Lys Val Trp Glu Glu Val Thr Pro Leu Thr Phe Ser Arg 

40 45 50 

Leu Tyr Glu Gly Glu Ala Asp He Met He Ser Phe 
55 60 



20 



35 



Ala Val Arg Glu His Gly Asp Phe Tyr Pro Phe Asp Gly 
65 70 75 



Pro Gly Asn Val Leu Ala His Ala Tyr Ala Pro Gly Pro 
25 80 85 90 

Gly He Asn Gly Asp Ala His Phe Asp Asp Asp Glu Gin 
95 100 

30 Trp Thr Lys Asp Thr Thr Gly Thr Asn Leu Phe Leu Val 

105 110 115 

Ala Ala His Glu He Gly His Ser Leu Gly Leu Phe 
120 125 



His Ser Ala Asn Thr Glu Ala Leu Met Tyr Pro Leu Tyr 
130 135 140 



His Ser Leu Thr Asp Leu Thr Arg Phe Arg Leu Ser Gin 
40 145 150 

Asp Asp He Asn Gly He Gin Ser Leu Tyr Gly Pro Pro 
155 160 165 

45 Pro Asp Ser Pro Glu Thr Pro 

170 



50 
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10 



15 



20 



25 



30 



35 



(2) INFORMATION FOR SEQ ID NO: 2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 83 amino acids 

(B) TYPE: amino acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: peptide 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 

Met Ala Thr Thr Pro lie He His Leu Lys Gly Asp Ala 
5 10 

Asn He Leu Leu Cys Leu Arg Tyr Arg Leu Ser Lys Tyr 
15 20 25 

Lys Gin Leu Tyr Glu Gin Val Ser Ser Thr Trp His Trp 
30 35 

Thr Cys Thr Asp Gly Lys His Lys Asn Ala He Val Thr 
40 45 50 

Leu Thr Tyr He Ser Thr Ser Gin Arg Asp Asp Phe Leu 
55 60 65 

Asn Thr val Lys He Pro Asn Thr Val Ser Val Ser Thr 
70 75 

Gly Tyr Met Thr He 
80 



(2) INFORMATION FOR SEQ ID NO: 3: 



(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 18 base pairs 
40 (B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



45 



50 



(ii) MOLECULE TYPE: DNA (genomic) 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 
GAAATGAAGA GTCTTCAA 18 
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10 



15 



25 



30 



35 



(2) INFORMATION FOR SEQ ID NO: 4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 18 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 4: 
GCGTCCCAGG TTCTGGAG 18 



(2) INFORMATION FOR SEQ ID NO: 5: 

20 (i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 28 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 



(ii) MOLECULE TYPE: DNA (genomic) 
(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 
ATACCATGGC CTATCCATTG GATGGAGC 28 

(2) INFORMATION FOR SEQ ID NO: 6: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 30 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
40 (D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 

45 

ATAGGATCCT TAGGTCTCAG GGGAGTCAGG 30 
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WHAT IS CLAIMED IS: 

1 . A process of screening compounds lo identify compounds that arc ligands 
that bind to a specific target molecule comprising the steps of: 

a) generating a first two-dimensional ^^N/^H NMR correlation 
spectrum of a ^^N-labeled target molecule; 
5 b) exposing the labeled target molecule to one or a mixture of chemical 

compounds; 

c) generating a second two-dimensional ^ ^N/ ^ H NMR correlation 

spectrum of the labeled target molecule that has been exposed to one 
or a mixture of compoimds in step (b); and 
10 d) comparing said first and second two-dimensional ^^N/^H NMR 

correlation spectra to determine differences between said first and 
said second spectra, the differences identifying the presence of one 
or more compounds that are ligands which have bound to the Urget 
• molecule. 

15 

2 . The process of claim 1 wherein the l^N-labeled target molecule is 
exposed to a mixture of chemical compounds in step (b)» further 
comprising the steps subsequent lo step d) of 

e) exposing the l^jsj-labeled target molecule individually to 
5 each compound of said mixture, 

0 generating a two-dimensional ^^N/ ^ H NMR correlation 
spectrum of the labeled target molecule that has been 
individually exposed to each compound; and 

g) comparing each spectrum generated in step f) to said first 
10 spectrum to determine differences in any of those compared 

spectra, the differences identifying the presence of a 
compound that is a ligand which has bound to the target 
molecule. 

3 . The process of claim 1 wherein the differences in the two-dimensional 
^^N/^H NMR correlation spectra are chemical shifts at particular ^^N- 
labeled sites in the target molecule and chemical shifts in protons attached to 
those ^^N-labelcd sites. 
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The process of claim 1 wherein the target molecule is a polypeptide. 

A process of determining the dissociation constant between a target 
molecule and a ligand that binds to that target molecule comprising the steps 
of: 

a) generating a first two-dimensional ^^N/ ^ H NMR correlation 
spectrum of a ^^N-Iabeled target molecule; 

b) exposing the labeled target molecule to various concentrations of a 
hgand; 

15 1 

c) generating a two-dimensional Nl H NMR correlation spectrum at 
each concentration of ligand from step (b); 

d) comparing each spectrum from step (c) both to the first spectrum 
from step (a) and to all other spectra from step (c) to quantify 
differences in those spectra as a function of changes in ligand 
concentration; and 

e) calculating the dissociation constant between the target molecule and 
the ligand from those differences according to the equation: 

Kn= ([P]o-x)([L]o-x) 

X 

where [Plo is the total molar concentration of target 
molecule; 

[L]o is the total molar concentration of ligand; and 
X is the molar concentration of the bound species detenmined 
according to the equation: 

A 

where 6obs and Sf^e are the chemical shift values for the 

target molecule determined at each concentration of 
ligand and for the target molecule in the absence of 
Ugand, respectively, and A is the difference between 
the chemical shift at saturating amounts of ligand 
and 6free. 
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6. The process of claim 5 wherein the target molecule is a polypeptide. 



7. 



The process of claim 5 further comprising the step of binding the labeled 
target molecule to a second ligand before step (a). 
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